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ABSTRACT 


Porou*  graphite  and  specially  prepared  aluminum  foams  (40  to 
80  percent  of  crystal  density)  were  investigated  in  an  attempt  to  deduce 
•  the  effect  of  such  material  parameters  as  particle  (or  pore)  shape, 

size,  and  size  distribution  on  response  of  the  materials  to  shock  loading. 
Results  suggest  that  the  strength  of  the  solid  matrix  and  the  porosity  are 
more  important  than  particle  geometry  per  se. 

Hugoniota  for  one  aluminum  and  one  graphite  foam  (both  about 
three-fourths  of  crystal  density)  were  measured  in  detail  up  to  about 
25  kbar.  In  the  pressure  and  porosity  range  studied,  the  P-V  Hugoniota 
of  the  foams  after  compaction  have  been  found  to  be  very  close  to  those 
of  the  solid  materials.  Recovered  shocked  specimens  of  aluminum  foam 
exhibited  densities  characteristic  of  solid  aluminum,  whereas  recovered 
shocked  specimens  of  ATJ  graphite  had  densities  very  close  to  their 
initial  densities. 

An  artificial  viscosity  digital  computer  code  was  adapted  to 
foams  and  shown  to  be  capable  of  predicting  approximate  shock  attenu¬ 
ation  behavior  of  porous  solids  as  measured  experimentally. 
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SECTION  I 


INTRODUCTION 


Research  on  foams  (distended  solids)  was  undertaken  at  Poulter 
Laboratories  in  19621  because  :t  was  recognized  that  such  highly  com¬ 
pressible  materials  should  be  effective  in  reducing  the  peak  pressure 
delivered  to  a  solid  boundary  by  a  given  impulse  at  the  opposite  surface 
of  the  foam  material.  Since  that  time  work  on  representatives  of 
several  classes  of  distended  materials,  including  plastics,  metals, 
ceramics,  and  graphite,  has  been  performed  at  Poulter  Laboratories 
under  U,  S.  Air  Porce  Weapons  Laboratory  sponsorship.1’3’4 

Although  all  the  foams  previously  studied  have  shown  similar 
behavior  in  a  gross  sense,  there  are  important  differences  in  their 
relative  effectiveness  as  countermeasure  materials  The  current  pro¬ 
gram  of  research  was  therefore  undertaken  to  study  two  types  of  foam 
(aluminum  and  graphite}  in  detail  in  an  attempt  to  (1)  discover  the  eifect 
of  varying  individual  material  parameters,  (2)  investigate  the  nature  of 
foam  Hugoniots,*  and  (3)  study  shock  attenuation  in  the  foems  as  a  basis 
for  comparison  with  theoretical  predictions. 


*  A*  used  in  this  report  the  term  '’Hugoniot"  is  applied  to  foams  to 
mean  the  locus  of  final  macroscopic  pressure-volume  or  pressure- 
particle  velocity  states  of  Jhocked  material  deduced  from  experi¬ 
mental  observations  and  invocation  of  the  so-called  "jump”  conditions. 
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SECTION  II 


THEORETICAL  BACKGROUND 

The  Hugoniot  equation  of  state  of  foams  ie  often  approximated  by 
the  elastic  locking-solid  model,  which  is  represented  in  the  pressure- 
volume  plane  in  Fig.  1.  For  this  model,  the  material  behaves  like  an 
isotropic  elastic  solid  below  a  certain  critical  pressure  P  .  Above  the 
critical  elastic  y/ave  amplitude  associated  with  the  point  (P^,  V^), 
pore  space  collapses  irreversibly,  and  the  material  locks  at  volume 
V  ,  so  that  it  behaves  like  an  incompressible  solid  by  comparison  with 
the  porous  material.  For  the  limiting  case  of  a  simple  locking  solid, 

Pe  *  0  and  Ve«Vo.  In  reality,  the  precise  shape  of  the  horizontal  por¬ 
tion  of  the  schematic  Hugoniot  of  Fig.  I  would  be  expected  to  depend 
upon  the  particular  foam  material  in  question,  but  should  ideally  be  a 
shallow,  rising  curve  with  positive  curvature. 

In  practice,  if  calculations  of  3hock  interaction  or  attenuation 
are  to  be  made,  the  vertical  portion  of  the  curve  is  often  replaced  with 
the  Hugoniot  of  the  solid  material  from  which  the  matrix  of  the  foam  was 
constructed.  Although  it  is  impossible  (on  thermodynamic  grounds)  for 
specimens  of  different  porosity  to  have  Hugoniots  v/hich  exactly  coincide 
in  a  region  of  the  pressure-volume  plane,  it  may  in  practice  be  a  very 
good  approximation  for  a  given  matrix  material  over  a  certain  range  of 
porosities  and  shock  pressures,  as  is  demonstrated  by  the  experimental 
results  discussed  in  this  report. 

Let  us  consider  a  shock  traveling  through  material  initially  in  a 
state  characterized  by  pressure  P^t  specific  volume  V  ,  and  total  in¬ 
ternal  energy  . 

Now  if 

E  =  E  (P,  V) 
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But  if  thermal  equilibrium  is  assumed  {as  an  approximation),  we  have 
(3E/ap)v  =  v/r.  where  T  is  the  so-called  Griineisen  coefficient;  and 


.  p 


(4) 


where  C  is  the  heat  capacity  at  constant  pressure  and  a  is  the  coeffi¬ 
cient  of  thermal  expansion.  Therefore 


dP 

dV 
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- EL  +  p  .  p 

vet  0 
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(5) 


Thus,  in  general,  the  slope  dP/dV  becomes  infinite  when  Vj/V=l  +  2/r 
and  becomes  positive  when  V^/V^  l  +  Z/T,  so  that  the  limiting  compres¬ 
sion  is  equal  to  1  +  Z/T  unless  the  material  has  an  unusual  equation  of 
slate  or  a  first»order  phase  transition  occurs.  In  other  words,  if  we  have' 
il  porous  solid  with  V0  >V#,  where  V8  is  the  initial  specific  volume  of  ,noa- 
porous  material,  then  the  limiting  compacted  volume  will  be  greater 
than  the  limiting  volume  for  initially  nonporous  material.  In  fact,  if  ■ 
V0/Vt>  1+  Z/T  the  Hugoniot  would  be  expected  to  have  positive  slope, 
since  a  small  pressure  P>Pfc  should  compact  the  foam  to  very  nearly 
solid  density  but  a  strong  shock  can  compact  it  no  more  than  to  a  volume 
of  V0/{l+2/r).  If  Vo/V#  =  1+2/r  we  have  in  principle  the  vertical 
Hugoniot  associated  with  the  ideal  locking  solid.  For  V /V^<l  +  2/r 
we  have  above  P#  the  more  familiar  Hugoniot  shape  commonly  associ¬ 
ated  with  the  shock  compression  of  initially  solid  materials.  These 
cases  are  illustrated  schematically  in  Fig.  2. 

For  purposes  'of  illustration  let  us  explore  a  specific  case.  *  A 
commonly  used  equation  of  8t..te  (Mie-Gruneisen)  for  solids  takes  the 
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whene  E^_  is  the  internal  energy  of  material  compressed  to  volume  V 
at  a  temperature  of  0°K  (i.e. ,  Ec  is  the  "compressional  part”  of  the 
total  internal  energy  £);  P  is  the  pressure  at  volume  V  and  tempera¬ 
ture  T;  and  A  is  an  undetermined  coefficient.  Now  if  for  simplicity 
the  heat  capacity  at  constant  volume  Cy  is  assumed  not  to  vary  with  T, 

then  £  =  £,.+  CVT,  and  we  have 

HR  E-E 

p  +  ^,r_£_  (?) 

wh«src  r  =  A/Cy  *  V(6P/3E)y  is  the  Gruneisen  coefficient. 

Using  the  conservation  of  energy  relation 
e-so  =  i(P0+P)(V-V) 
with  E^  and  set  arbitrarily  to  zero  yields 

3E  »v  /  ,  y  N 

V~sr-*T\  *  -v-r\**r-v-J  (8> 

along  the  Hugonlot,  so  that 

&£ 


where  Of  =  (5V/t»T)p/V  is  the  coefficient  of  thermal  expansion  and 
K  *  (3VfCiP)j./V  is  the  bulk  isothermal  compressibility.  At  room  tem¬ 
perature  and  atmospheric  prcisure,  this  yields  F*-  2.  1  for  aluminum. 


For  solids,  r  will  generally  fall  between  1  end  4.  Thus,  the  critical 
distention  for  "cross-over"  of  the  Hugoniot  should  be  somewhere  be¬ 
tween  V  /V.  *1.5  and  V  /V.  *  3.  As  the  pressure  increases,  V  /V 

d  I  0  •  O 

approaches  the  limiting  value  1+2/T.  Therefore,  the  limiting  wlue  of 
V  1c  proportional  to  the  initial  specific  volume  V0 ,  provided  that  T 
remains  constant  (as  it  should  for  the  pressure  range  considered  in  this 
report).  In  actuality,  T  may  be  expected  to  decrease  with  decreasing 
volume  and  with  increasing  temperature,  but  should  still  be  on  the  order 
of  1  to  2  even  at  very  high  pressure. 

The  theoretical  basis  for  foams  serving  as  effective  pressure 
attenuators  is  covered  in  detail  in  Reference  1  and  will  be  only  quali¬ 
tatively  summarised  hero.  Foams  are  effective  as  countermeasure 
materials  because,  first  of  ail,  a  Hugoniot  which  resembles  that  of 
Fig.  1  implies  a  Urge  ratio  cf  rarefaction  velocity  to  shock  velocity. 

A  shock  front  travels  at  a  veloelr  t  proportional  to  the  square  root  of  the  - 
slope  of  tha  straight  line  connect!)  g  the  initial  and  the  shocked  states  in 
the  P-V  plane.  The  sum  of  the  particle  velocity  and  longitudinal  sound 
velocity  in  the  compressed  material  behind  the  shock  is  proportional  to 
the  square  roof  of  the  elope  of  the  adiabat  through  the  shocked  state. 

The  first  means  of  attenuation  may  be  thought  of  aa  caused  by  rarefactions 
overtaking  the  shock  puUe  from  the  rear.  Therefore,  foams  are  espe¬ 
cially  effective  because  they  provide  a  Urge  rarefaction-to-oompactioa 
velocity  ratio  (see  Fig.  1),  Alternative  ways  of  looking  at  the  same 
effect  in  terms  of  energy  dissipation  or  of  "momentum  spreading"  lead 
to  the  tame  conclusions. 

A  very  important  advantage  can  be  giined  by  the  presence  of  a 
forerunner  wave.  For  shock  pressures  between  certain  limits  (i.e„, 
between  P%  and  Pc  In  Fig.  1),  the  shock  front  splits  up  into  two  waves, 
the  first  of  theee  waves  having  amplitude  P  .  Since  momentum  and 
energy  must  be  conserved,  attenuation  of  tlx  second  (slower)  wav*  will 
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j  be  enhanced  as  the  amount  of  material  included  between  the  two  wave- 

/ 

fronts  increases.  The  rate  of  attenuation  will  also  be  proportional  to 
•  the  momentum  per  unit  volume  of  this  included  material.  Thus,  it  is 

desirable  to  have  the  velocity  difference  between  the  waves  as  high  as 
possible  and  to  have  the  value  of  Pe  as  large  as  possible,  provided  that 
the  critical  damage  pressure  in  the  structure  to  be  protected  will  not  be 
exceeded  when  the  forerunner  wave  reaches  the  structure -foam  interface. 


i 


/ 


.  ! 
f 

f 


SECTION  IU 

SHOCK  ATTENUATION  CALCULATIONS 

1.  Description  of  Method 

Calculation*  of  shock  propagation  and  interaction  may  be  made 
precisely  using  the  method  of  characteristics;  however,  this  method 
soon  becomes  unwieldy,  when  used  for  any  but  the  simplest  problems  and 
is  difficult  to  organize  into  a  digital  computer  code.  As  an  alternative, 
the  artificial  viscosity  method  (Q-method)  of  von  Neumann  and  Richtmyer* 
is  readily  adapted  to  computer  calculations  but  has  the  possible  disad¬ 
vantage  of  spreading  the  shock  front  ovor  a  finite  distance.  The  Q-method 
uses  finite  difference  equations  ih  real  space  and  time  to  represent  dif¬ 
ferential  equations  describing  the  flow?*  7  The  artificial  viscosity  Q 
is  added  to  the  pressure  in  the  difference  equations  to  remove  all  dis¬ 
continuities  in  the  flow. 

The  Q-method  has  been  used  for  all  shock  attenuation  calculations 
reported  here. 

The  value  of  Q  used  was  given  by 


«„Cip]/v 


(11) 


where  and  ConR  are  arbitrary  constants  which  were  set  equal  to 

Z'O  and  0. 1,  respectively;  C  is  the  Local  acoustic  velocity,  V  is  the 

■P 

local  specific  volume,  and  &u  is  the  particle  velocity  increment  (i.e. , 
local  particle  speed  difference  between  adjacent  cells  of  constant  mass). 
It  should  be  noted  that  Q  as  given  above  has  a  negligible  vaU.c  outside 
the  shock  region. 

The  computations  were  carried  out  on  a  Burroughs  B5500  digital 
comp'  from  an  ALGOL  program  developed  by  John  O.  Erkman  and 
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adapted  to  porouc  solid*  by  George  Duvall.  Th'j  essential  steps  in  the 
calculation  are  summarized  in  Appendix  I,  where  the  entire  program 
has  also  been  reproduced. 

2.  Equations  of  State 

The  equation  of  state  information  has  been  programmed  as 


follows: 


2  lb/ft*  Styrofoam 


Th#  relations 


[V  -V  (V-V  ){V-Va)  *| 

vrb~^J 

/•v  <  V  <  V 

{ •  *■ 

^0  <P<P 


■where  P4  “l.  debar,  *  31.25  cm*/g,  V,  *  0.9«2  ere*/*,  and 

b  *  0. 1:  and. 


P  *  A 


VG  -V 
—  -  1  * 


■(H' 


-  «  •  » 

/  * . 

Ip‘*  r  '  .  »;• 


:  r  V*  V 


wteere  A  *  90.5  kbart  5M  200  kbsir,  and  0.  9524  aa*/g  .{based 

en  the  data  of  Wagner  et  *1-*  for  1.05  g/o?1  polystyrene)  have  been  . 
"'••d  in  the  STYROFOAM  procedure  faee  Appendix  2). 

6061-T&  Aluminum  Flyer  Plate 


The  relationa  . 


»"  U£'*%**.***&*tr9i>»<  t 


P  =  1044 


.P  =  794 


for  0  s  P  £  6.  4  kbar 


+  1.6  for  6.4  <P*  31  kbar 


have  been  given  by  Lundergan  and  Herrmann*  and  were  adapted  to  the 
elastoplastic  model  using  the  ELASTGPLASTIC  procedure  developed  by 
Erkman10  (see  Appendix  I). 

Quartz 

•The  relation 


P  =  869.5  (1-V/V  }  kbar,  (15) 

0 

where  V  *  3.774  ca’/j,  has  been  devised  from  the  data  of  Graham 
o 

et  fcl.  and  used  in  the  WATER  procedure  (see  Appendix  I)  to  represent 
the  equation  of  state  of  quartz  under  both  shock  and  rarefaction.  The 
value  0.5728  cm/psec  has  bc»-  used  for  G  of  Eq.  (11). 

Porous  Targets 

For  the  porous  aluminum  and  graphite  targets  a  new  procedure 
called  PHASE  TRANSITION  which  allows  hysteresis  tc  exist  in  a 
compression-rarefaction  cycle  nas  been  written  by  George  Duvail.  The 
available  P-V  information' for  the  parous  targets  and  the  -t3sumption3 
made  hero  are  best  described  with  reference  to  Fig.  3. 

The  experimental  results  reported  in  Section  VI  give  values  for 

the  amplitude  P  of  the  forerunner  wave  in  the  target  material  and  for 

the  corresponding  specific  volume  V^.  These,  with  the  initial  volume 

V  ,  give  'he  end  points  of  a  quadratic  which  has  negative  slope  and 
o 

positive  curvature  for  V  s  V  *  ^ ,  provided  that  0<bs  1 : 
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Along  AB 


-(V-V  1  (V-V  )(V-V  /  ' 

f  (V)  =  P.  ~rr~  +  B - S - ! 

1  1  Yv0  iv.-vf 

L  x  o  -i 


where  V  e  0.478  cm1/*,  V  =  0.474  cm%,  P  -  0.446  kbar  and 
o  l  °  1 

b  *  0.  1  for  the  2. 1  g/cm  aluminum  foam  studied  and  V  =  0.  580  cm3/g, 

o 

V  »  0.  576  cm3/g,  =  0.  27  kb  ir,  b  =  0.  1  for  the  ATJ  graphite. 
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FK).  3  ASSUMED  P-V  DIAGRAM  IN  PHASETRANSITIGN  PROCEDURE 


The  quadratic  form  for  the  P-V  curve  in  this  region  was  chosen 
for  the  following  reason:  If  a  wavefront  is  stable,  the  loco*  of  state*  ia 
the  P-V  plane  followed  by  a  mass  element  passing  through  ths  wavefront 
is  «  straight  line,  l.e. ,  the  effective  normal  stress  P^*P+Q»  a(V-V0). 
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The  gradient  in  the  wavefront  then  adjust*  itself  so  that  Q  *  *(V  -  V„)  -  P. 
If  the  curvature  of  the  P,-V  relation  is  aero,  the  corresponding  steady- 
state  value  of  Q  is  aero;  however,  the  rise  time  of  the  wave  will  depend 
upon  the  value  of  Q,  as  discussed  in  Section  VI.  In  order  to  allow  some 
freedom  in  adjusting  the  rise  times  to  match  experiment  (Section  VI-4), 
b  has  been  given  a  small  positive  value.  II  rise  times  are  of  no  interest, 
the  value  of  b  can  be  set  equal  to  aero. 

The  curve  BCD  connecting  the  elastic  curve  AB  with  the  solid 
curve  GDE  was  chosen  in  the  Murnaghan  form  (but  with  arbitrary  K): 


f,M-  v-f[(4)  •■] 


Values  of  parameters  used  were  Px  =  0.46  kbar,  Vt  =.0.474  cm*/g, 

K  *  12.06  and  R  *  2.0  for  2. 1  g/cm*  aluminum  foam,  and  Px  =  0.23  kbar, 
Vx  *  0.576  cm*/g,  K  *  2.909,  and  R  *  12.0  for  ATJ  graphite  foam. 
Equation  (17)  is  a  curve  with  positive  curvature  passing  through  Px,  Vx 
which  allows  partial  compaction  between  (P4,  Vx)  and  (Pf,  V#)  and  insure* 
compaction  to  crystal  density  at  (Pt,  V#)(see  Section  VI),  where  the  re¬ 
sulting  valuta  of  P  and  V.  are  4.  5  kbar  and  0. 367  cm*/g  for  the  aluminum 
foam,  and  7,5  kbar  and  0.433  cm*/g  for  die  graphite  foam.  The  values 
of  Px  and  Vx  were  determined  experimentally  (Section  VI),  while  the 
values  of  the  other  constants  were  adjusted  to  be  consistent  with  die 
attenuation  experiment*  of  Section  VJ-4  (rather  than  directly  determined), 
and  therefore  are  highly  subject  to  the  effects  of  any  errors  in  the  value 
of  Px «  ........ 

The  equation  used  for  curve  GBE  was  of  the  form 


where  A  *  729  kbar,  B  *  1717  kbar,  and  Vq  *  0. 369  cm*.  ^  for  the 
2.  1  g/cm*  aluminum  foam,  and  A  *  193  kbar,  B  *  2855  kbar,  and 


i\.. 


I 


*  **•  f  r'/  <. 


*-  «*\  ”  4 

*>  •,*»  !( 
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Vq  »  0.444  cmJ/|  for  the  ATJ  graphite.  The  values  of  the  constants 
for  aluminum  were  inferred  on  the  basi3  of  extrapolation  of  the  measure¬ 
ment!  of  Al'tshuler  et  al. 11  and  normal!  nation  of  che  data  of  Lundergan 
and  Herrmann*  on  6061-T6  aluminum.  In  the  caoe  of  ATJ  graphite,  the 
constants  were  chosen  by  normalizing  the  data  of  Doran. x* 

The  compression  process  is  assumed  to  correspond  to  a  Hugoniot 
represented  by -the  curve  ABODE  specified  by  Eqs.  (16)  to  (18).  If  ex¬ 
pansion  occurs  before  point  B  is  reached  in  compression,  it  is  assumed 
to  occur  along  BA.  If  expension  occurs  after  point  D  is  reached,  it  Is  • 
assumed  to  occur  along  £DG. 

For  lack  of  better  information  it  is  assumed  that  if  expansion  . 
occurs  from  any  point  C  along  BCD,  it  occurs  along  a  path  CF  which  is  - 
chosen  ac  cor  dir.;  tn  the  following  prescription; 


Along  CF: 


Y-V 

P  s  f  (v)  “  P  1  +  - — 

■  int'  f  c  V  -V 
c  co 


(V-V  )(V-V  )  I 
+  b - £ - — 

(V  -V  )*  | 

c  co  J 


where  (V^,  P^)  is  the  point  on  BCD  where  compression  stopped  and 
expansion  began,  V  is  obtained  from  the  simple  proportion: 


V  -V 
CO  Q 

V  -V 
X  o 


v  _v 
c  ; 

V  -V 
a  '  i 


and  the  other  quantities  have  the  same  values  as  specified  previously. 
Thie  recipe  is  simple  and  easy  to  compute  and  it  farces  the  expansion  ■ 
curve  to  coincide  with  AB  when  C  coincides  with-B.  It  effectively  re¬ 
places  Eq.  (18)  below  point  D,  but  the  difference  between  the  two  curves 
is  small  between  D  and  G.  It  ij  furtaer  assumed  that  if  a  reversal  from 
expansion  to  compression  occurs  while  the  state  point  lies  on  FC,  com¬ 
pression  will  proceed  to  point  C  and  then  along  CDE  as  before. 


I 

I 


for  both  compression  and  rarefaction  wave*.  The  equations  used  along 
AB>  BD,  and  DE  were  the  same  as  the  compression  curves  discussed 
above,  with  the  exception  that  along  AB  the  value  b  =  0. 1  was  used  for  . 
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compression  but  fe  *  0  was  used  for  rarefaction,  so  as  to  avoid  '  ! 

artificially  introducing  a  rarefaction  shock.  (A  flow  chart  is  given  in  | 


Fig.  25,  Appendix  I.)  "  \ 
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SECTION  IV  . 

material  development...!.. 

2.  Specimen' Fabrication 

To  perform  reliable  experiment*,  the  target  materials  must 
have  reproducible  behavior  and  must  be  homogeneous  on  a  scale  small 
compared  with  their  thickness.  In  the  case  of  carbon  and  graphite 
materials  we  have  of  necessity  relied  upon  the  commercial  supplier. 

ATJ  graphite  was  chosen  for  most  of  the-  work  because  it  provides  the 
best  combination  of  desired  properties  (i.s. ,  homogeneity  and  yield 
strength).  Other  types  of  graphite  arid  carbon  were  used  for  the  param¬ 
eter  variation  studies;  although  they  were  less  desirable  materials, -- 
X-rity  transmission  ar.d  volumetric  density  measurements  indicated 
them  to  be  reasonably  homogeneous,  with  the  exception  of  GA  Coke- 
carbon.  All  graphite  or  carbon  foams  were  essentially  of  the  "open  cell" 
type  (i.  e. ,  they  were  not  impervious  to  the  passage  of  gases). 

In  the  case  of  porous  aluminum  specimens,  a  considerable  amount 
of  developmental  work  was  done.  Suitable  aluminum  foams  are  not 
commercially  available.  The  ~nly  commercial  aluminum  foam  used  was 
MO-AK(see  Appendix  II),  which  was  rather  inhomogeneous,  although  its  -  - 
low  density  and  high  yield  strength  indicate  it  to  be  very  promising  as  a 
shock  pressure  attenuator  (see  Section  VI). 

The  major  problem  in  making  aluminum  foams  stems  from  the  - 
herd,  nearly  impervious  oxide  layer  which  forms  on  all  aluminum  - 
surfaces.  We  have  employed  three  techniques  for  producing  homogeneous 
foams  in  the  nominal  density  range  of  40  to  80  percent  of  crystal  density 
(depending  on  particle  shape  and  pressure  of  pressing).  These  techni¬ 
ques  are;  (1)  hot  pressing  of  aluminum  powders;  (2)  cold  pressing  of 
aluminum  powder*  with  subsequert  sintering;  and' (3)  repeated  hot  •  V- 


pressing  and  sintering  of  a  mixture  of  aluminum  powders  and  thin- 
walled  silica  microballoons. 

The  hot  pressing  was  accomplished  by  heating  (to  450-570° F) 
aluminum  powders  of  known  particle  size  distribution  in  a  hardened 
steal  die  and  then  pressing  to  the  desired  density  (using  700  to  Z0,  000 
psi). 

For  technique  2,  the  aluminum  powder  was  cold  pressed  (9500 
to  13.000  psi)  to  essentially  the  desired  final  density.  Thin  specimens 
ware  pressed  into  solid  aluminum  rings  (with  1/8-inch  walls)  for  ease 
of  handling.  The  sample  disks  were  placed  in  a  stainless  steel  box  on 
graphite  pads.  The  box  wits  put  in  a  furnace  and  an  argon  flow  was 
maintained  through  the  box.  The  temperature  of  the  furnace  was  raised 
from  room  temperature  to  1000°  F  ovnr  a  period  of  five  hours  and  held 
one  hour  at  1000°  F;  power  to  the  furnace  was  then  shut  off.  The  furnace 
cooled  to  200»300°F  in  about  five  hours,  during  which  time  the  argon 
flow  continued.  ,  t 

In. the  case  of  aluminum  flake,  the  powder  did  not  cold  bond.  It 
wae  therefore  necessary  to  add  a  small  amount  cf  wax  to  make  compacts 
that  would  hold  together  well  enough  to  be  handled.  These  compacts 
were  heated  more  slowly  to  give  the  wax  a  chance  to  fume  off.  Maximum 
sintering  temperature  waa^l000°F; 

Under  technique  3,  aluminum  powder  and  silica  microballoons 
(walls  about  2  microns  thick)  were  mixed  in  an  elbow  blender  and 
the  mixture  was  placed  in  a  graphite  pressing  die  in  a  pot  furnace.  The 
temperature  was  raised  to  900° F  over  a  period  of  five  hours.  The  die 
was  thee  removed  and  the  mixture  was  pressed  at  about  500  psi.  The 
material  was  then  returned  to  the  furnace  ar.d  heated  to  1000° F,  held 
on*  hour,  removed  and  pressed  O  500  psi,  returned  to  the  furnace  at 
1100°F  for  one  hour,  pressed  agtin  to  590  psi,  and  then  ajlowejl  to  . 
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cocl  to  room  temperature.  Densities  achieved  indicated  that  about 
’ialf  the  porosity  was  due  to  the  microballoons  and  half  to  ordinary  gas- 
filled  pores. 

Except  for  the  commercially  obtained  MD-AK  all  foams  were 
essentially  of  the  "open  cell"  type. 

With  Technique  2  it  was  easier  and  less  expensive  to  control 
specimen  quality  than  with  Technique  1,  although  there  is  no  inherent 
superiority  of  the  method.  Specimens  were  made  in  the  shape  of 
3 Jf*  inch-diameter  ~dtsks  of  various  thicknesses  and  were  subsequently  ‘ 
machined  to  else  if  necessary  (e.g. ,  in  the  case  of  the  quasi-static 
compression  samples).  Technique  3  was  used  only  for  foams  of  density 
*60  percent  of  crystal  density. 

2.  Metallography  and  Quality  Control  • 

The  metallographic  method  involved  saw-cutting  specimens' 
{oriented  parallel  io  and  perpendicular  to  the  pressing  direction)  from 
the  sintered  aluminum  foam  samples.  Special  care  was  taken  to  avoid 
deformation  of  the  foam  in  the  regions  to  be  examined.  The  samples 
were  vacuum-impregnated  with  epoxy  resin  and  catalyst  (Armstrong 
C-7).  Surfaces  were  prepared  by  wet  sanding  on  silicon  carbide  coated 
papers,  100  through  600  menh,  hand  lapping  with  30-micron  alumina; 
and  final  polishing  with  6-micron  diamond.  In  some  cases  it  was  neces¬ 
sary  to  reimpregnate  the  polished  surface  with  epoxy  and  repeat  the  . 
sending  and  polishing  operations.  The  loose  powder  samples  were 
mounted  either  in  Ducite  or  in  epoxy,  depending  on  the  nature  of  the 
particles.  The  particle  sixes  observed  in  the  fcam  samples  correspond 
to  the  screen  fraction  used  and  to  the  particle  size  and  shape  observed, 
in  the  powder  mounts;  hence  it  is  reasonable  to  tssume  that  the  struc¬ 
ture  of  the  foam  observed  is  real  and  not  the  result  of  gross  distortion, 
cold  flow,  etc.  .  ‘ 
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Types  of  powders  used  and  the  resulting  foams  are  ohovsn  in  the 
photomicrographs  of  Figs.  5  and  6.  No  differences  were  noticed  be¬ 
tween  photomicrographs  taken  at  orientations  parallel  to  and  perpendic¬ 
ular  to  the  direction  of  pressing.  The  small  dark,  spots  in  the  particles 
of  Figs.  6a  and  6c  are  probably  polishing  pits. 

.  In  addition  to  the  metallographic  examination,  selected  specimens 
were  checked  for  uniformity  of  density  by  either  X-ray  transmission  or 
sectioning  and  weighing.  Longitudinal  acoustic  velocities  were  measured 
(see  Sections  V  andVl)  for  each  specimen  as  a  means  of  quality  control. 

The  aluminum  powder  used  in  the  foam  for  which  the  detailed 
Hugoniot  measurements  to  25  kbar  were  made  (see  Fig.  6a,  b,  c)  and  ou 
which  the  attenuation  studies  were  performed  was  the  "chunk"  powder 
type  MD801-  (see  Appendix  II)  of  crystal  density  2.71  g/cm*.  The  manu¬ 
facturer's  spectroscopic  analysis  disclosed  the  composition  to  bet 

A1  99.5%  ■ 

•  •  Fe  0,15  .  ' 

Si  0.06 

...  Zn  0.02  ■  ..v  •*  ... 

.  .  .  .  O  .  0. 27 
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SECTION  V 


EXPERIMENTAL  TECHNIQUES 

1.  Shock  Wave  Measurements 

a.  Technique  for  Measurement  of  Hugoniot  Point* 

All  Hugociot  measurements  (including  the  parameter  vari¬ 
ation  studies)  employed  the  impedance-match  method  (see  Fig.  7).  The 
projectile  (made  from  2024-T351  aluminum,  a  material  with  a  reasonably 
well  known  Hugoniot*)  was  accelerated  to  a  velocity  uA  and  impacted 


FKk  7  ILLUSTSiTION  OF  IMPEDANCE-MATCH  METHOD  FOR 
MEASURE  TARGET  HUGONiOT 


*  Extrapolation  of  the  data  of  Fowleg**  for  2024-T4  aluminum,  in  conjunc¬ 
tion  with  other  data  at  f,  13,  and  17  kbar. 


the  target  material.  The  initial  state  in  the  projectile  is  represented 
by  point  A  of  Fig.  7,  and  the  initial  state  in  the  target  material  is  rep¬ 
resented  by  point  O.  The  conservation  relations  (jump  conditions)  re¬ 
tire  that  upon  impact  both  materials  must  reach  the  pressure-particle 
velocity  state  represented  by  point  B.  Conservation  of  momentum  re¬ 
quires  that  across  a  single  shock 


B 


P 

o 


B 


*oU. 


'where  p  Is  the  Initial  density  of  the  target  material  and  is  the 
velocity  of  a  shock  of  amplitude  characterised  by  P_,  and  u  in  the 
material.  Thus,  measuring  the  shock  velocity  experimentally  and 
knowing  the  initial  density  of  the  target,  a  line  through  point  O  can  be 
drawn  with  slope  p  U  .  At  point  B  this  line  intersects  the  known 
projectile  Hugoniot  drawn  from  f  .int  A.  Therefore,  point  B  must  lie 
on  the  Hugoniot  of  the  target.  Using  different  projectile  velocities  (i.e. 
various  value*  of  u^),  we  can  trace  out  the  entire  Hugoniot  for  the  tar¬ 
get  material. 


If  Pg  of  Fig.  7  is  less  than  Pc  of  Fig.  1,  a  single  shock  is  un¬ 
stable  and  the  pressure  pulse  in  the  target  will  consist  of  two  fronts, 
oti«  of  amplitude  Pjj  (P,  of  Fig.  1),  and  the  other  of  amplitude  Pg.  The 
determination  of  the  target  Hugoniot  is  thon  more  complicated,  but  if 
Pj.«  Pg  no  great  error  will  be  introduced  by  adopting  the  following 
procedurej  as  illustrated  in  Fig.  8. 

1.  Measure  1*^  with  a  quart*  gtge,  as  described  below. 

2.  Measure  (T'  ^  and  a  draw  line  through  O  with  slope 

p  (U  )  (point  E  it  as  yet  undetermined), 
o  *  oe 

3.  On  the  basis  of  other  cross-curve  data,  obtained  at  higher 
pressures,  estimate  the  cross  curvu  from  the  line  found  in 
step  2  to  point  Q  measured  in  step  1.  This  fixes  point  E. 


Su0P^*Vu»5«  particle  velocity  tt.w 

PIG,  I  ILLUSTRATION  OP  PROCEDURE  FOR  HUGONIOT  MEASUREMENT 
WHEN  A  DOUBLE  WAVE  STRUCTURE  EXISTS- 


4.  Calculate  Pj.  from  the  jump  condition*  and  draw  line  EB. 

5,  Proceed  u  in  single  shock  css®. 

If  p  ••  p  ft,  uncertainty  in  the  estimated  cross  curve  became* 

2  » 

very  Important’*  iteration  may  help.  A  serious  problem  is  also  intro¬ 
duced  because  (U  will  be  ia  error  due  to  reverberationo  of  the  fore- 
runner  between  the  quarts  ga  ye  and  the  main  w.rvefront.  If  Pg  »  P£ » 
however,  the  uncertainty  in  :ha  cross  curve  is  of  rather  minor 
importance,  and  not  oniy  will  the  reverberation i  of  the  forerunner 

.  /  -  zr 
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occur  over  a  relatively  email  percentage  of  the  total  transit  time  o-er 
which  (Us5gg  ie  meaeured,  but  each  reverberation  results  in  only  a 

•  mall  perturbation  in  (U  }__  . 

8 

Implicit  in  this  entire  development  if  the  assumption  that  steady 
state  has  been  reached.  Furthermore,  since  shock  velocities  are 
measured  as  averages  over  the  transit  time  in  the  specimen  (see  below) 
considerable  error  may  be  introduced  if  steady  state  has  not  been 
achieved  in  a  time  short  compared  with  the  total  transit  time.  The  situ¬ 
ation  is  clarified,  however,  by  repeating  experiments  with  different 
target  thicknesses  but  the  same  projectile  velocity.  In  the  experiments 
reported  herein,  forerunner  pressures  (P^)  have  not  fejen  found  to  be  as 
reproducible  as  desired;  however,  confirmation  of  Hugoniot  points  for 
Pg  >>  P#  ha*  been  excellent  {see  Section  VI). 

Except  where  noted,  all  experiments  employed  quarts  gages  as 
the  primary  sensing  element.  The  quarts  gage  is  an  electroded  (with 
guard  ring)  X-cut  synthetic  quarts  disk.  Current  from  X-cut  quarts 
crystals  may  be  used  as  a  sensitive  indicator  of  stress -time  profiles, 
at  discussed  in  detail  by  Graham  et  al. 11  The  current  amplitude  and 
it*  time  dependence  are  functions  of  the  dielectric,  pieaoelectric,  and 
mechanical  properties  of  shcck-loaded  quarts'. 

The  normal  stress  difference  between  die  two  faces  of  the  disk 
is  given  initially  fay  the  relation  .*• 


RAKV 


where  V  is  the  potential  drop  across  the  tofcil  load  resistance  R 
through  which  the  pieaoelectric  current  flows;  .1  is  the  thickness  of  the 
disk;  A  la  the  active  electrode  area,  K  is  an  exp-irimenully  determined 
current  coefficient,  and  is  the  shock  velocity  in  quarts.  Values  of 
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U  and  K  are  given  by  Graham  et  al. x  at  a  iunction  of  stress.*  Rather 
than  using  their  suggested  average  values  for  each  pressure  range,  we 
have  chosen  to  draw  best-fit  curves  to  their  data;  the  discrepancy  is 


small,  however.  For  the  quartz  gages  usod  in  our  experiments^ 

A  *  0.201  in.  *,  A  •  0.376  in. ,  and  R  =  51. 0  0,  so  that  Fig.  9  may  be 


used  to  determine  the  normal  pressure  at  the  specimen-quartz  inter¬ 
face  from  measurement  of  the  voltage  across  R. 


as  a  function  of  time.  As  the  shock  propagates  through  the  quartz,  the 

curve  of  Fig.  9  will  shift  continuously  in  the  direction  of  higher  voltages 

*  The  values  of  K  above  IJ1  kbar  and  of  above  25  kbar  show  some 
scatter,  however,  so  that  computed  pressures  arc  probably  less 
reliable  at  the  higher  pressures. 
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for  the  ume  pressure,  perhaps  reaching  as  much  as  a  7-perceut  in¬ 
crease  for  transit  of  a  20-kbar  shock,  so  that  pressure-time  profiles 
may  be  distorted  and  corrections,  should  be  made  to  avoid  errors  in 
interpretation.  Since  the  magnitude  of  the  shift  depends  on  shock  am¬ 
plitude  within  the  quartz  disk,  results  for  the  second  wave  of  a  double 
wave  profile  should  not  be  affected  much  by  the  passage  of  low-amplitude 
forerunner  waves.  As  described  below,  the  gage  was  surrounded  with 
potting  material  which  approximates  the  shock  impedance  of  quartz,  so 
that  for  the  shock  attenuation  experiments  the  useful  lifetime  of  the 
gage  could  be  extended  by  reconstructing  approximate  profiles  (i.e. , 
since  the  voltage  is  proportional  to  AP,  the  value  of  P  at  the  specimen- 
quarts  interface  could  be  estimated  by  adding  the  stress  value  found  one 
transit  time  earlier — that  is,  1. 66  fisec  earlier). 

Average  shock  velocities  were  determined  from  the  time  inter¬ 
val  between  tilt  pin  closure  (see  below)  and  the  time  at  which  the  pro¬ 
file  for  each  shock  reached  half  its  total  associated  pressure  increment 
(see  diagram  accorr  anying  Table  VI,  Appendix  HI).  The  use  of  this 
"half -wave  point"  for  measuring  velocity  can  be  argued  on  the  basis  of 
the  profile  shapes  (see  Fig.  21)  and  has  the  fortunate  advantage  of  pro- 
vidiuF  the  least. opportunity  for  random  reading  errors.’ 

b.  Experimental  Details 

Target  Preparation.  The  prime  constituents  of  the  target 
assembly  are  the  specimen,  the  adapter  ring,  the  tilt  pine,  the  sensing 

element  (quarts  gage),  and  the  radial  velocity  pins  (see  Figs.  10  and  11).  j 

Each  specimen  was  lapped  to  a  tolerance  of  0.02  mm  to  0.0002  mm 

I 

(depending  on  specimen  material  and  projectile  velocity)  and  then  j 

attached  to  a  prefabricated  molded  support,  which  was  thick  enough  to  ’  { 

prevent  bowing  of  the  specimen  (see  Fig.  10).  This  support  was  con-  i 

structed  from  a  mixture  of  Armstrong  C-?  epoxv  and  glass  beads  j 


A  " 
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(40  H  to  75  H  size  rang''  in  a  proportion  adjusted  to  approximate  the 
shock  impedance  of  the  quart*  gage.  A  thin  layer  of  quick-setting 
epoxy  was  applied  to  the  lapped  surfaces  of  the  molded  support  and 
the  specimen,  which  were  held  in  contact  by  a  uniform  pressure  until 
the  cure  cycle  for  the  adhesive  was  complete. 

The  four  tilt  pin  clearance  holes  were  drilled  at  90°  intervals 
using  a  guided  No.  51  drill;  the  tilt  pin  wires  (copper.  No.  22)  were 
bent  and  inserted  through  the  rear  face  of  the  specimen  until  they  ex¬ 
tended  about  i  mm  beyoru.  the  impact  surface.  The  bared  portion  of 
each  tilt  pin  was  supported  near  the  center  of  the  clearance  hole  by  the 
insulation  (which  covered  the  conductor  for  approximately  half  its 
length  In  the  specimen).  Quick-setting  epoxy  was  used  to  fill  to  excess 
the  remaining  volume  between  the  pins  and  the  walls  of  the  clearance 
holes,  so  as  to  give  good  bonding  for  positive  location  and  mechanical 
support.  On  the  rear  surface  of  the  specimen  the  conductors  were  also 
potted  with  quiek- setting  epoxy  for  support. 

After  the  pins  had  been  permanently  located  with  respect  to  the 
specimen,  the  excess  potting  material  and  pin  lengths  were  removed  "by 
filing  and  lapping  until  the  leading  surface  of  -each  pin  lay  within  0.  002 
mm  of  the  plane  defined  by  the  impact  surface  of  the  specimen  as  deter¬ 
mined  by  a  dial  gage.  Any  small  deviation  from  this  plane  was  noted 
and  considered  when  the  shot  data  were  reduced. 

The  quarts  gage  was  clamped  firmly  in  place  on  the  rear  surface 
of  the  specimen,  and  a  small  amou  of  viscous  quick-setting  epoxy 
was  applied  around  its  circumference.  The  resulting  bond  held  the 
gage  in  place  during  subsequent  potting  and  formed  a  seal  that  did  not 
allow  any  potting  material  to  seep  between  the  gag*  and  the  specimen. 

To  check  that  the  gage  was  in  good  contact  with  the  specimen,  a  com¬ 
bined  thickness  measurement  was  maclu  across  the  specimen  end  gage 
and  was  compared  with  the  sum  of  their  individual  thicknesses. 


After  the  attachment  of  an  electrical  ground  and  the  necessary 
signal  cables  and  bias  resistor,  the  assembly  was  potted  with  the  pre¬ 
viously  described  epoxy-glass  bead  mixture.  When  the  potting  com¬ 
pound  had  cured,  electrical  checks  were  made  to  determine  if  thexe 
were  any.  shorted  pins  or  faulty  connections.  The  assembly  was  placed 
in  a  carefully  mavhined  recess  in  the  adapter  ring  and  fixed  in  place 
with  Shell  Ep>.i  (No.  91  IF)  adhesive,  which  provided  a  vacuum-tight 
seal. 

The  adapter  ring  was  also  used  to  locate  the  three  retractable 
radial  velocity  pins.  These  pins  were  located  by  using  high-precision 
guide  holes  in  the  adapter  ring  to  define  a  plane  parallel  to  that  of  the 
specimen  at  an  accurately  measured  distance  of  about  6  mm.  The  dis-' 
charge  of  these  pins  upon  contact  with  the  projectile  was  detected  and 
recorded  as  described  in  detail  below. 

Projectiles.  The  ?.  6-cm-long,  hollow  2024-T35 1  aluminum 
projectiles  used  in  all  Hugoniot  shots  (see  Fig.  10)  were  fabricated  on 
a  lathe  and  subsequently  placed  in  a  close-fitting  guide  .ring  for  machine 
lapping.  Final  squareness  and  flatness  of  the  projectile  head  were 
generally  held  to  within  0.002  mm.  Projectile  heads  were  either  6  mm 
or  i2  mm  thick.  For  the  shock  attenuation  shots  it  was  desired  to  im¬ 
pact  the  target  with  a  short-duration,  high-amplitude  pressure  pulse 
rather  than  the  flat,  long-duration  pulse t  provided  by  thick  projectile 
head*.-  To  this  end,  0. 0160-inch-thick  6061-T6  aluminum  flyer  plates 
ware  prepared  In  the  form  of  2.  3fc-inch-diameter  disks.  Each  flyer  - 
plate  was  mounted  on  a  1 -inch-thick  piece  of  low  density  (2  lb/cu  ft) 
styrofoam,  which  in  turn  was  mounted  on  the  head  of  a  modified  2024 
aluminum  projectila  (see  Fig.  12).  Care  Ws  taken  to  mount  the  flyer 
plat*  perpendicular  to  the  projectile  walls. 
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FIG,  12  PROJECTILE  ASSEMBLY  FOR  SHOCK  ATTENUATION  STUDIES 


By  u«e  of  this  flyer  plate'  assembly,  it  was  possible  to  achieve 
relatively  plane  impact  with  the  desired  pulse  length  {~0.  15  jisec)  and 
input  momentum  densities  (~0.  5-Q.75  x  104  dyne-sec/cm*).  The 
styrofoam  support  could  produce  only  a  very  minor  perturbation  on  the 
shock  profile,*  and  the  2024  aluminum  projectile  could  play  no  role  - 
whatever  during  the  times  at  which  profile  measurement*  were  made. 

All  profile  measurements  employed  the  quartz  gage  techniques  described 
below.  "* 

Target  Alignment  Procedure.  It  is  extremely  important  to  • 
minimise  the  relative  tilt  1v tween  projectile  and  target  beeause  shock 
wave  velocities  are  typically  an  order  of  magnitude  greater  than  the 
associated  particle  velocities  even  for  foams  at  the  shock  strengths  used 
in  thla  study. 

*  The  effect  of  the  styrofnun  was  implicitly  taken  into  account  in  the 
attenuation  calculations  (see.  Section  HI);  a  comparison  of  calculated 
shock  profiles  with  and  vithout  the  styrofoam  nupport  is  giyon  in 
Section  VI. 
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.The  target  assembly  was  mounted  on  the  muzzle  end  of  the  gun, 
the  position  of  the  target  being  controlled  by  four  close-fitting  dowel 
pins, and  three  fine-threaded  alignment  screws.  A  reference  plane  per¬ 
pendicular  to  the  axis  of  the  gun  barrel  was  provided  by  the  alignment 
tool  (see  Fig.  11).  The  screws  were  used  to  rotate  the  free  surface  of 
the  specimen  into  this  plane  without  breaking  the  vacuum  seal. 

The  alignment  tool  is  a  46-cm  brats  cylinder  that -closely  fits 
th*  •■ore  of  the  gun  barrel.  The  precise  clearance  required  is  controlled 
by  using  three  adhesive  cellophane  tape  strips  spaced  120°  apart  on  the 
circumference  and  extending  the  length  of  the  cylinder.  A  short  cello¬ 
phane  flap  extends  onto  the  ends  of  the  tool.  The  build-up  of  these  strips 
involves  the  use  of  r wo  or  three  layers  of  tape,  each  ir.  turn  being  cut 
to  length  from  material  5  mm  wider  than  the  previous  layer. 

Four  adjustable  hardened  steel  pine  and  a  light  spring  tension 
ground  wire  are  mounted  on  a  Micarta  head  attached,  with  an' epoxy  ad¬ 
hesive  to  one  and  of  the  alignment  tool.  The  steel  pine  are  adjusted  and 
h«od-lapped;  to  within  0.002  mm  of  each  other.  Electrical  conductors 
extend  from  the  pins  out  through  the  breach  plug  of  the  gun  to  a  12-volt 
power  supply  and  four  indicator  lamps.  Contact  of  each  pis  with  ths 
(conductive)  specimen  completes  a  circuit  that  turns  on  a  specific  indi¬ 
cator  lamp.  Contact  pressure  is  supplied  by  a  weak  spring  located  be- 

'  .  » 
twesn  the  breach  plug  and  a  long  rod  "handle"  behind  the  brass  cylinder. 

The  target  alignment  procedure  (see  Fig.  11)  consisted  of  (1) 
sliding  the  alignment  tool  down  the  gun  barrel  until  it  contacted  the 
target,  (2)  evacuating  the  system,  (3)  turning  the  thrte  fine-threaded 
alignment  screws  until  contact  between  the  adapter  riag  and  the  vacuum 
manifold  was  broken  (the  ends  of  th?  screws  and  the  "0"-ring  then  pro-’ 
vided  total  support  for  the  target  assembly),  and  (4)  adjusting  the  align¬ 
ment  screws  until  all  four  indicator  lamps  were  glowing  uniformly. 
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The  alignment  wag  checked  by  venting  the  gun  barrel  to  atmospheric 
pressure  and  rotating  the  alignment  tool  by  120°,  240e ,  and  360°  (re¬ 
evacuating  each  time). 

If  all  four  lamps  were  lit  for  all  three  positions  of  the  alignment 
tool,  the  target  was  pe-pendicular  to  the  axis  of  the  gun  barrel  to  withih 
about  0.002  mm.  The  failure  of  all  four  lamps  to  light  simultaneously 
after  rotation  of  the  alignment  tool  was  an  indication  that  the  target  vat 
not  fla»  to  the  required  tolerances.  The  deviation  from  flatness  could  be 
estimated  by  the  amount  that  the  alignment  screws  must  be  turned  in 
order  to  light  the  unlit  lamps  in  the  rotated  configuration. 

The  tilt  pin  arrangement  described  above  allowed  the  measure¬ 
ment  of  both  magnitude  and  orientation  of  tilt.  Each  pin  was  charged  to 
+  22v,  and  discharged  upon  contact  with  the  (conductive)  projectile.  To 
display  the  occurrence  times  of  electrical  pin  closures  unambiguously, 
a  binary  coded  system  was  employed  (see  Fig.  13).  The  tilt  pins  were 
charged  by  a  lfiF  capacitor  that  formed  part  of  an  RLC  circuit.  When 
oach  pin  was  grounded  it  discharged  through  a  different  resistor,  re¬ 
sulting  in  a  different  current  and  ultimately  a  different  measured  voltage 
across  the  50  0  output  resistor.  *  The  signals  from  the  four  plus  were 
conveyed  through  a  single  cable  and  displayed  on  an  internally  triggered, 
cathode -modulated  (see  below),  fast-rise  oscilloscope  oweeping  at  0. 1 
ftsec/cm  (see  Fig.  14a).  As  tho  flatness  of  impact  is  improved,  the 
voltage  steps  become  more  and  more  coincident  until  only  one  large 
step  is  seen  (Fig.  14b). 

Projectile  Velocity  Measurement.  The  first  Hugoniot  shots  fired 
(Shots  11,  117  to  11,316)  employed  a  linear  array  of  electrically  charged 

*  The  resulting  signals  are  given  dative  amplitudes  of  1,  2,  4,  and  8 
(see  record  of  Fig.  14a)  by  use  •  •{  the  circuit  resistors  specified  in 
Fig.  13. 


FIG.  13  BINARY  CODER  CIRCUIT  FCR  DISTINGUISHING  ELECTRICAL 
PIN  CLOSURES 


contact  pins  which  extended  through  the  side  of  the  gun  barrel  and  were 
discharged  through  the  projectile  (see  Fig.  10)  for  deducing  projectile 
velocities.  The  measured  times  at  which  these  pins  discharge  were 
used  to  compute  projectile  velocity  and  acceleration;  the  results  were 
than  extrapolated  to  determine  the  velocity  at  time  of  impact.  Subse-  • 
quent  Hugoniot  shots  employed  the  more  accurate  radial  velocity  pin 
syotem  pictured  In  Fig.  10.  The  radial  velocity  pins  are  retractable 
brass  rods  with  0. 25-mm-diameter  "whisker"  ends  that  are  electrically 
charged  to  +  lOOv  and  operate  through  a  circuit  similar  to  that  of 

13.  They  lie  in  a  plane  6  rr.m  frc  n  the  impact  surface  of  the  tar¬ 
get  and  may  be  extended  radially  into  the  path  of  the  projectile  after  the 
target  has  been  aligned  (see  Figs.  10  and  11).  The  final  position  of  each 
pin  Is  controlled  by  a  carefully  machined  shoulder  in  each  guide  hole 
(see  Fig.  10).  Timing  signals  are  prode,  ed  when  the  charged  pins  ar,'  ' 
grounded  upon  physical  contact  with  the  c inductive  projectile  (which, 
in  turn,  is  coupled  to  ground  through  the  jan  barrel).  One  of  the  radls.i 
pins  is  0.05  mm  farther  from  the  target  and  connecieo  to  ground. 

The  radial  pin  signals  were  displtyed  on  an  intsrnally  triggered, 
cahtode -modulated,  fnst-rlse  oscilloscope  sweeping  at  0.  1  (S sec/cm. 
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FIG.  14  OSCILLOSCOPE  RECORDS  SHOWING  TILT  PIN  CLOSURES 
AND  CATHODE  MODULATION  TIME  MARKS 
(a)  Torgtt  no!  pr«ci«»ly  aligntd;  ord*r  of  pin  closure  it  2,  8,  1,  4. 

.(b)  Torgtt  oligrtd;  pin  cloiurtt  crt  •  sssntially  timultontous. 

This  allowed  small,  positioning  errors  in  the  pins  to  be  averaged  out 
and  also  allowed  tilt  o£  the  projectile  to  *be  taken  into  account  in  the 
velocity  computation*. 

The  first  radial  pin  closure  signal  also  started  a  10-mc  counter 
and  was  displayed  on  a  raster  oscilloscope.  The  first  tilt  pin  closure 
stopped  the  counter  and  was  displayed  oi  the  raster  oscilloscope.  The 
information  thus  obtained,  when  used  in  conjunction  with  the  cathode- 
vnodulated,  binary-coded  oscilloscope  records,  was  sufficient  to 
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calculate  to  within  about  0.  01  pse c  the  time  interval  between  the  pro¬ 
jectile  arrival  at  the  plane  of  the  radial  pint  and  the  impact  with  the  tar¬ 
get  (see  below). 

Time  Correlation!.  It  it  important  to  establish  a  well  defined 
time  bate,  to  that  correlations  can  be  made  between  events  recorded 
by  different  instruments.  If  an  otcillotcope  trace  has  a  fiducial  mark 
occurring  every  0.  SO  jjsoc,  then  two  separate  events,  such  as  the  clos¬ 
ing  of  a  tilt  pin  and  the  occurrence  of  a  quarts  gage  signal,  may  be  in¬ 
dividually  recorded  on  different  oscilloscopes  cet  to  trigger  internally 
or  through  uo me  delay.  These  oscilloscopes  may  be  set  to  sweep  at 
0.  1  (iitc/cm  to  that  much  detail  is  resolved.  If  the  relative  occurrence 
times  of  events  are  known  to  the  nearest  half-microsecond  {perhaps  from 
a  slow -sweeping  otcillotcope),  the  fiducial  marks  may  be  used  to  pro¬ 
vide  the  necessary  accurate  correlation,  provided  the  fiducial  marks 
themselves  are  accurately  controlled. 

The  cathode  modulation  timing  system  described  below  minimize* 
errors  due  to  variations  in  oscilloscope-  sweep  speeds  and  trace  linearity, 
and  Is  capable  of  providing  time  correlations  between  oscilloscopes  to 
within  0.01  /jses.  Care  must  be  taken.,  however,  to  employ  the  norma, 
precautions  to  prevent  high-frequency  noise  from  appearing  in  the  data 
signal  traces  (1.  e.,  use  shielded  cables,  ground  oscilloscope  chassis 
etc  ). 

The  time  marks  originate  from  an  E-H  Model  120D  pulse  ganar- 
aior*  capable  of  supplying  pulses  0.  01-jisec  wide.  -The  pulse  generator 
waa  controlled  by  an  external  high-stability,  temperature-compensated, 

2-  nc  crystal  oscillator.  The  timing  pulses  were  used  to  modulate  the 
cathode  bias  voltage  of  the  oscilloscope  CkT's  (aee  Fig,  15}.  Pulse 

*  Manufactured  by  E-H  Research  Labor  .tories,  Ins  .  Oakland, 
California. 
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K3TE.WI)  ALL  CASCX3  ARE  3ifl  COAX 

(2)  CALIBRATED  Slfl  TEBM‘NATtOh 
RESISTORS  ARE  PROVIDED  AT 
POINTS  A,B,G,0,E,ANOF 


PIN  Dl$<'  LAT 

FIG.  15  BLOCK  DIAGRAM  OF  INSTRUMENTATION  FOR  RECORDING 
AND  CALIBRATING  DATA  SIGNALS 

amplitudes  were  adjusted  so  that  the  time  mark*  were  represented  by 
blanking  of  the  oscilloscope  traces  at  0.  5-psec  Intervals.  To  allow  cor-- 
relation  of  different  signals,  the  electrical  cables  labeled  in  Fig.  15 
we  chosen  to  be  of  the  same  type  and  have  lengths  such  that 

7B  *  JC  -  BCf  a  JID  -  BCD  *  JKJT  -  BCDE,  (22) 

The  lengths  of  the  other  cables  shown  in  Fig.  15  were  not  important, 
since  the  calibration  procedure  compensates  for  their  effects. 

The  calibration  procedure  involved  connecting  all  the  dotted 
cables  as  indicated  in  Fig.  15  and  removing  the  termination  resistors  at 
A,  B,  C,  and  D,  while  leaving  uie  oscilloscope  sveep  speed  settings  and 
all  other  cables  and  delay  units  the  same  as  durii the  gun  shot.  Oscillo¬ 
scope  records  made  with  the  system  in  this  configuration  appear  as  in 
Fig.  16 a,  where  a  cathode  modulation  timing  mark  occurs  at  time  tft 
on  the  trace  ,  ar.d  the  maximum  of  the  blanking  pulie  which  produced  the 
time  mark  is  displayed  at  time  t  .  The  interval  (t.  -l^}  ropresents  the 
sum  of  the  interna!  delays  in  the  oscilloscope  and  tie  time  delaya 
associated  with  external  delay  unit*  and  the 
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extra  length  of  cable  through  which  the  pulse  must  travel  between  the 
cathode  and  the  vertical  input  terminal  of  the  oscilloscope.  A  compari¬ 
son  of  the  value  of  (tj  -  t0)  for  each  oscilloscope  was  used  to  adjust  the 
shot  data  to  a  common  time  base.  *  Since  the  time  marks  were  uaed  at 
0.  5-fisec  intervals  only,  a  10-mc  sine  wave  was  used  to  check  linearity 
of  the  oscilloscope  sweep  between  time  marks.  A  typical  oscilloscope 
record  from  a  quartz  gage  shot  using  cathode  modulation  time  marks  is 
shown  in  Fig.  16b.  Shots  11,269  to  11,  799  employed  the  time  mark 
system  described  above. 

2 .  Measurement  of  Longitudinal  Sound  Velocities 

Figure  17  shows  a  block  diagram  of  the  sound  velocity  measuring 
system.  A  G.6-fi*ec  square-wave  pulse  from  a  pulse  generator  is  con¬ 
verted  into  small-amplitude  approximately  sinusoidal  mechanical  waves 
of  1  me /sec  frequency  by  a  *«t  of  two  Independent  piezoelectric  (P2^T-4) 
crystal*  contained  in  the  mercury  delay  line  and  the  aample  holdsr.  The 
mechanical  waves  traverse  their  respective  media  and  are  reconverted 
lmo  electrical  aignala  by  the  aecond  set  of  PZT  crystals.  These  ''output" 
signals  may  be  viewed  on  an  oscilloscope;  their  difference  is  found  by 
the  Type  1A1  amplifier.  The  length  of  the  Hg  column  is  carefully  adjusted 
so  that  this  difference  equals  *ero  (i.e. ,  Signal  A  ~  Signal  B);  the  system 
is  then  nulled  and  the  time  delsy  through  tire  mercury  and  the  sample 
are  the  same- 

,  The  longitudinal  velocity  of  sound  in  mercury  (1.446  nun/f wee  at 
256G)  Is  accurately  known  and  quite  insensitive  (-0.  000465  c/°C) 

to  reasonable  temperature  changes;  thus  it  makea  an  ideal  medium  for 
comparison  purposes.  The  "null  method"  for  determining  the  velocity 
of  sound  in  a  solid  material  requires  that: 

*  If  the  cable  lengths  satisfy  Eq.  (22),  a  common  time  base  is  estab¬ 
lished  by  simply  Mbtraetliv,  the  appropriate  value  of  (t,  *  te)  from 
the  times  indlcab  i  by  the  cathode  modulation  time  marks  in  esch 
shot  record. 

43 


/ 


PULSE  I 
GENERATOR) 


WBffl 

H«  oelay  line 


2 

SAMPLE 


SCOPE 

TRIGGER 


AMPLIFIER 

TYPE  IAI  PlUfl- IN-UNIT 


OA-MOtSSI-lM 


where 


FIG.  17  BLOCK  DIAGRAM  OF  SYSTEM  FOR  MEASUREMENT 
OF  LONGITUDINAL  SOUND  VELOCITY 


t  »  delay  time  through  the  mercury 

Hg 

t  a  delay  time  through  the  sample 


where 


dH  «  path  length  through  the  mercury 
v,,  •’  1.44$  mm/jMec  at  25ttC. 


t  ..  »  time  de.\y  at  aero  thickness  of  the  mercury 

oHg 

column  {constant)  -  <  ' 

d^  “  path  length  through  the  specimen  - 

t-  «  time  delay  at  sero  thickness  of  the  specimen 
0  8  ■* 

holder  (constant) 

v  *  the  longitudinal  (sound)  velocity  being  measured. 
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SECTION  VI 


EXPERIMENTAL.  RESULTS  AND  DISCUSSION 

1.  Parameter  Variation 

The  parameter  variations  are  illustrated  in  part  by  the  photo¬ 
micrographs  shown  in  Section  IV.  A  wide  range  of  densities  and  par¬ 
ticle  {or  pore)  shapes,  sites,  and  size  distributions  was  investigated. 
Unfortunately,  ther  e  is  too  much  scatter  in  the  data  (see  Appendix  ID) 
to  allow  specific  conclusions  as  to  the  effect  of  any  single  parameter 
alone  in  the  above  list;  however,  certain  general  observations  can  be 
made. 

a.  MD-AK  aluminum,  which  has  the  strongest  matrix  (solidi¬ 
fied  molten  aluminum,  rather  than  the  weaker  sintered  bonds  between 
particles--see  Figs.  6  and  ?)  has  the  largest  precursor  amplitude  and 
speed,  despite  its  low  density  (1.45  g/cms).  Tht  situation  is  not  en¬ 
tirely  clearcut  because  MD-AK  has  large  spherical-type  pores  rather 
than  small  irregular  ones  as  do  the  sintered  materials.  However,  our 

preliminary  results  indicate  the  pore  shape,  size,  and  size  distribution 

» 

to  be  of  less  importance  than  matrix  strength. 

b.  For  both  graphite  and  sintered  aluminum  foams,  the  pre¬ 
cursor  amplitude  and  speed  decrease  with  increasing  distention,  but 
experiments  on  aluminum  foams  show  that  distention  alone  may  not  be  as 
important  as  matrix  strength,  which  may  be  expected  to  depend  upon 

the  manufacturing  process;  • 

’  c.  In  the  pressurs  and  porosity  range  studied  for  each  type 
of  foam,  dynamic  data  show  that  the  "compacted"  volumes  for  pres¬ 
sures  above  a  few  kilcfcars  are  essentially  those  of  the  solid  materials, 
regardless  of  pjrticle  (or  pore)  shape,  size,  or  si»e  distribution, 
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resulting  in  the  plots  of  Fig*.  18  to  20  (see  Section  VI-C).  Possible  ex¬ 
ception*  are  MD-AK  (for  which  we  do  not  have  reliable  Hugoniot  point*), 
aluminum  flake  fane  *hot  only;  «ee  also  Footnote  a  to  Table  III),  and 
GA  Coke  carbon(one  ahot  only;  *ee  also  Footnote  c  to  Table  IH).  Points 
representing  compacted  ATj  graphite  show  no  dependence  on  orientation 
and  have  not  been  labeled  with  orientation  in  the  figures.  (Orientations 
are  given  in  Appendix  13. ) 

d.  Dynamic  elastic  limits  scale  approximately  with  quasi¬ 
static  elastic  limits  and  lie  generally  within  a  factor  of  three  of  the 
quasi-static  values.  Some  scatter  is  present  in  the  data,  however. 

a.  Precursor  velocities  are  generally  higher  in  foams  which 
hav#  higher  elastic  limit*.  Quasi-static  yield  street  mey  prove  to  be 
a  reasonable  indicator  of  relative  precursor  velocities  for  a  given  type 
of  material. 

f.  Measured  longitudinal  acoustic  velocities  (see  Section  V-B) 

7 

tend  to  be  qualitative  indicator*  of  elastic  wave  velocities,  but  are  gener¬ 
ally  larger  than  the  elastic  wave  velocity,  as  might  be  expected  for  foams 
(sea  Tables  V  to  VII  in  Appendix  HI). 

It  should  be  noted  that  Lcckyer  has  investigated  nondestructive 
methods  for  evaluation  of  graphite  materials  (including  ATJ  graphite) 
and  found  correlations  between  such  quantities  as  longitudinal  sound 
velocity,  bulk  density,  electrical  resistivity,  tensile  elastic  modulus, 
ultimate  tensile  strength,  and  dynamic  elastic  modules  in  flexure.  Thus, 
it  is  reasonable  to  suppose  that  some  inexpensive,  nondestructive  tests 
mey  prove  useful  indicators  of  shock  forerunner  amplitude*  r.xd  speeds 

I.  Terminal  Observations  on  Recovered  Materiel 

All  of  the  shots  fired  were  designed  for  the  purpose  of  naking 
dynamic  measurements  on  the  foams.  However,  in  many  case?  pieces 
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of  the  target*  which  had  experienced  Loth  the  ir.  .,al  shock  and  the  shock 
reflected  from  the  foam-quartz  interface  were  recovered  in  a  condition 
suitable  for  measurements  of  density. 


Results  are  shown  in  Table  I,  where  the  column  headed  "P  " 
contains  the  range  of  maximum  pressures  of  inflected  shocks  experi¬ 
enced  by  material  for  which  the  reported  recovered  dc  cities  were  actu¬ 
ally  measured.  In  the  case  of  "chunk"  aluminum  and  1  TJ  graphite,  ma 
values  reported  for  recovered  density  are  based  on  the  *vera;x  for 
several  recovered  targets.  In  the  case  of  all  other  materials  the  value 
reported  is  based  upon  data  for  only  one  or  two  shots.  Table  II  contains 
recovered  density  information  for  quaai-statically  compressed  {with 
•  Ides  confined)  specimens  to  serve  as  a  basis  for  comparison.  These 
data  should  be  considered  is  only  approximate  typical  values.  In  many 
cases  the  data  are  based  on  only  one  experiment. 

The  results  are  very  Striking.  Sintered  aluminum  foams  evi¬ 
dently  "lock"  at  close  to  solid  density,  whereas  graphite  'cams  {partic¬ 
ularly  ATJ)  recover  their  porous  character.  Both  kinds  of  foam  are  of 
the  open  cell  type  (ae  determined  by  a  vacuum  integrity  test),  eo  that 
the  behavioV  of  graphite  as  contrasted  with  that  of  aluminum  ie  rather 
interesting.  The  MD-AK  aluminum,  on  the  other  hand,  is  a  "closed 
cell"  type  of  foam,  so  it  weald  be  expected  that  gases  would  be  trapped 
in  the  pores  during  shock  compression.  It  is  not  known,  however,- 
whether  the  lower  density  of  the  recovered  MD-AK  material  is  due  to 
"spring-beek"  or  whether  the  foam  was  never  fully  compacted-' 

Those  results  for  graphite  and  sintered  aluminum  foams  have 
been  used  in  estimating  release  curves  for  the  shock  attenuation  calcu¬ 
lations  (Section  III);  results  are  given  in  Section  VI-D.  The  "spring- 
back"  phenomenon  in  graphite  might  be  of  further  interest  because  it 
offers  tha  possibility  of  a  countermeasure  material  which  can  function 
for  more  than  one  "Wt.  '* 

49 


Table  I 

DENSITY  DATA  FDR  RECOVERED  SHOCKED  ?OAM3 


Material 

Estimated* 

t>  (kbar) 

■ex 

Initial 

Density 

(t/cm*) 

Typical 

Recovered 

Denalty 

(g/ca1) 

Alunlnua: 

Chunk 

(146-420  n) 

IS,  -..,60 

2.09 

2.68 

Milled 

(260-640  n) 

8-10 

1.06 

2.6 

Spherical 
(8-106  >i) 

8-10 

1.74 

2.0 

Spherical 
(44-106  m)  ♦ 
Mlc/oballoona 
(44-106  ») 

8-10 

1.63 

2.6 

MD-AJC 

6-10 

1.45 

2.3 

MO -AS 

80-70 

1.48 

2.8 

Graphite: 

ATJ 

10,  .,.,80 

1.73 

1.7 

Grade  60 

<v  28 

1.02 

1.68 

Eatlaated  at  quarts-fOM  Interface  aa  described  la  Section  71-3 


Table  II 

RESIDUAL  DENSITIES  OF  QUASI -STATICALLY  COMPRESSED  FOAMS 


Material 

P  (kbar) 
max'  7 

Initial 

Density 

(r/c»») 

— 

Typical 

Density 

W 

(s/c*») 

Typical 
Recovered 
Density  in  Die 
(«/«») 

Alurainu*; 

Chunk  ' 

(149-420  p) 

1.6 

2.09 

2.5 

2.5 

Milled 

(260-140  p) 

1.6 

1.03, 

1.40 

2.5 

2.45 

Spherical 
(6-106  u) 

l.S 

2.11 

2.56 

2.56 

Spherical 
(44-106  jl)  .+ 
Mlcroballeons 
(44-106  * 

1.6 

1.26 

2.38 

2.33 

MD-AK 

4.3 

1.46 

2.42 

2.31 

Graphite : 

ATJ 

1.6 

1.72 

2.12 

1.73 

Grade  60 

1.0 

1.02 

1.94 

1.59 
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Hu ?or.i ot»  of  Foams* 


The  plot*  of  Fig*.  18  to  20  sre  comprised  of  point.*  representing 
macroscopic  pressure-particle  velocity  states  for  various  graphite  and 
sintered  aluminum  foams.  In  all  cases  the  Hugcniot  points  were  deter¬ 
mined  by  the  impedance -match  method  as  discussed  in  Section  V,  and 
only  the  most  reliable  points  available  have  been  plotted.  Typical  voltage- 
time  profiles  measured  with  quarts  are  shown  in  Fig.  21.  Foam  targets 
of  substantially  different  thicknesses  yielded  the  same  Hugoniot  (i.e. , 
yielded  essentially  the  same  shock  velocity  for  the  same  impact  condi¬ 
tions),  indicating  that  the  shock  velocities  used  for  determining  the  Hugo- 
niots  must  have  been  vary  close  to  the  steady-state  velocltits.  Thick¬ 
headed  projectiles  were  used  for  the  Hugoniot  shot*  and  the  ‘‘flat-topped" 
profiles  measured  in  quarts  testify  to  the  fact  that  the  peak  pressure  at 
the  shock  front  was  not  being  attenuated  before  transit  through  the 
target.  Rarefactions  from  the  sides  of  the  target  would  also  have  changed 
the  ‘'flat-topped"  appearance  and  would  have  been  visible  in  the  guard- 
ring  record a  first. 


As  a  basic  for  interpretation,  smooth  curves  based  on  theoretical 
models  and  totally  Independent  of  the  measured  data  are  drawn -on  the 
same  axes  in  Figs.  If  to  20.  These  smooth  curves  ars  derived  from 
the  relations  of  conservation  of  maso  and  momentum  across  a  shock  front 
(jump  conditions)  and  taka  the  form, 


p  ■  p  ♦ 

i 


V  -V 
l  p 


where  P(kbar),  Up(mm/|isec).  and  V^(cm*/g)  are  die  pressure,  par¬ 
ticle  velocity,  and  specific  volume  associated  with  the  state  behind  the 
shock*  and  P^  (kbar),  (mw/jiMe),  and  V  (cm8/g)  are  the  pressure, 
psrticle  velocity,  and  specific  volume  associated  with  the  state  ahead 


*  Sea  loo  trot*  o»  f.  t. 


MATICU  VCLOCITV  — 

ICLS-VCLOCITY  DIAGRAM  FOR  2.1  «/«.*  "CHUNK" 
.  f-moe»h  cwrv««  end  i points  alonp  quo/tz  Hvpontot 
Utcrib*^  in  Hxt. 


of.  the  shock  (i.e.,  either  the  initial  "rest"  state  or  the  state  behind  the 

forerunner  wave,  as  the  case  may  be).  Curves  labeled  "compacted  to 

compressible  solid"  were  drawn  using  values  of  P,  u  ,  and  V  obtained 

P  P 

from  the  game  data  used  in  computing  Eq.  (16)  for  solid  graphite  and 


for  solid  aluminum,  with  allowance  made  for  any  forerunner  wave.  For 
purposes  of  comparison  the  curves  labeled  "compacted  to  incompres¬ 
sible  solid"  we~e  drawn  assuming  the  foam  to  collapse  to  the  atmos- 
•  pheric  pressure  crystal  density  and  to  remain  at  this  density  regardless 
oC  pressure  (rigid  locking-solid  model  of  Fig.  1).  Smooth  curves  are 


merely  labeled  "compactcd"in'cases  where  both  types  are  indistinguish¬ 
able  for  the  scales  employed.  Values  of  P^t  u^,  and  crystal  density- 
used  in  obtaining  the  curves  are  shown  in  Table  III. 


FIG.  19  PRESSURE,  PARTICLE-VELOCITY  DIAGRAM  FOR 

VARIOUS  ALUMINUM  FOAMS.  Sawoth  cuvet  on d  mmber*f 
point!  elect]  \JOr1i  Hvgonit*  weco  eeopvtoi  at  described  in  tcurt. 
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FIG.  30  PRESSURE  PARTICLE-VeLOCITr  DIAGRAM  FOR  GRAPHITE 

AND  CAR#ON  FCAMf .  Smooth  «t”ve*  end  number #«  point*  along  quart* 
Hugoolet  war*  campy 1*4  a*  AoteriM  In  text. 


(blShDTf.O  11,533  DOUBLE-WAVE 
PPC^'-E  'N  2  09  g/cms  ALUMINUM 
FOA’V 


Id)  SHfil  II.22C  DOUBLE-WAVE 
PPO^l.E  If.  ATj  GRAPHITE 


'fl  S“CTf.C  II. US  DOUBLE -WAVE 
PFCF  .E  Ifli  50>.ID  2024  - T35I 
aluminum 


.  1 16  FROM  QUAR“2 

tei*en-i,»ir:i  imerfc'<  os  described  in  Section  V-l 
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Tabls  III 

V ABASES  OF  QUANTITIES  USED  IN  DRAWING  .SMOOTH 
CURVES  OF  FIGURES  18  TO  20 


Material 

Initial  Density 

Forerunner 

Crystal  Density *• 

(g/e**) 

Pt(kbar) 

Ujfmm/nsec) 

(«/c*3) 

Aluminum 

Chunk 

2.09 

0.41 

0.015 

2.71 

Flake 

1.77 

0.0 

0.0 

2.71* 

Spherical 
+  micro- 
balloons 

1.33 

0.1 

0.004 

3.71 

(ND-AK)b 

1.45 

l.o 

c.oa 

2.71 

Spherical 
&  milled 

1.05 

0.0 

0..0 

2.71 

Graphite 

• 

■ 

ATJ 

1.72 

0.33 

0.03 

2.20 

Grade  CO 

1.02 

0.1 

0.01 

2.25c 

2.20 

AG  SR 

1.0« 

0.2 

0.01  ' 

2.25c 

2.30 

CKJ 

0.94 

0.0 

0.0 

2.20 

Carhop 

OA  Coke 

1.61 

0.0 

0.0 

CPA 

1.00 

0.1. 

■  ■  -  .1 

o.°o;i 

l.*- 

The  value  2.7  b/cb*  Bay  be  slightly  large  if  the  wax  waa  not  co»- 
plataly  -removed  during  heat  treatment  of  the  flak*  powders  (see 
Section  JY). 

^Th#  values  for  MD-AK  are  Included  bore  for  purposes  of  coops risen 
although  the  experimental  data  (see  Table  Til,  Appendix  tii)  did 
not  allow  precise  onough  measurement  of  Hugoniwt  points  %o  re¬ 
quire  the  ipclusior.  of  t  ,x>th  curves  for  MD-AK  in  the  figures. 

Cfhe  value  J.Xi  n/cv-k  1*  Jcr  natural  graphite,**  whereas 

the  valus  S,3Q  s/oa3  tr.  «  typical,  initial  dcnrlty  for  pyrolytic 
graphite,1* 

Shla  v»1m  lit  atmr'1  *  ^sonatile  density.  Possible  values 
range  trrm  1-  /cm*. 


5? 


Reference  to  Figs.  18  and  20  shows  that  in  the  cases  where 
several  reliable  data  are  available  the  "compacted  to  compressible 
•olid"  curves  are  essentially  best-fit  curves  to  the  experimental  data. 
This  result  indicates  that  the  pressure-volume  states  of  the  foams 
after  compaction  lie  indetectably  close  to  those  of  the  solid  materials. 
Data  points  below  a'b'ut  2  or  3  kbar  are  unreliable  for  some  materials 
because  of  uncertainties  associated  with  the  forerunner  wave  (see 
Section  V-A),  but  have  been  included  in  the  figures. 

These  results  are  reasonable  for  the  distentions  and  pressures 
involved,  in  light  of  the  discussion  cf  Section  U.  There  is  come  un¬ 
certainty  introduced  into  the  smooth  curves  by  the  extrapolation  or 
normalisation  of  data,  and  particularly  in  the  choice  of  the  appropriate 
crystal  density  for  the  graphite  specimens,  but  the  uncertainties  are 
not  large;  the  effect  of  changing  the  assumed  density  of  solid  graphite 
is  illustrated  in  Fig.  20.  The  results  obtained  lend  confidence  to  the 
assumed  P-V  Hugoniots  cf  Section  III,  with  regard  to  the  crystal  density 
curves  ODE  (Figs.  3  and  4;. 

A  very  important  point  is  raised  by  the  indicated  cross  curves  in 
Figs.  18  to  23.  Numbered  points  along,  the  quart*  Hugoniot  correspond 
to  r  assures  indicated  by  quart*  for  shots  in  which  the  foam  Hugoniot 
point  of  the  same  number  was  measured.  The  primed  points  represent 

c 

expected  pressures  in  quart*  for  the  given  foam  Hugoniot  point  on  the 
assumption  thst  compacted  foam  behaves  like  the  initially  solid  material. 
When  total  compaction  has  occurred  the  compacted  foam  should  behave 
like  the  solid  materia),  since  the  results  discussed  above  show  thst  the 
P-V  states  of  compacted  foam  end  initially  solid  material  lie  sc  close 
together  The  indications  from  the  quart*  voltages  that  the  cross 
curves  are  considerably  more  shallow  implies  he  dubious  result  of  s 
compressibility  for  the  "fully  compacted"  foam  much  greater  than  that 


for  the  initially  solid  material,  i.e.,  high  retained  porosities  (generally 
five  to  ten  percent)  would  have  to  be  presenl  even  at  the  higher  pressures, 
in  contradiction  with  the  Hugoniot  measurements;  furthermore,  this 
porosity  would  have  to  be  retained  after  one  strong  shock  but  removed 
after  two  weaker  ones. 

In  order  to  clarify  the  situation,  shots  were  fired  using  a  1/4-in.  - 
thick  2024-T351  aluminum  plate  between  the  foam  and  the  quartz;  others 
were  fired  using  Manganin  gages*7' 14  in  this  configuration.  The  results 
of  these  shots  showed  the  original  quartz-inferred  pressures  to  be  in 
error  and  the  assumed  cross  curves  (primed  points)  to  be  correct. 
Further,  shots  with  quartz  gages  on  solid  2Q24-T351  aluminum  targets 
indicated  pressures  in  quartz  that  were  in  accord  with  the  impact  condi¬ 
tions  and  known  Kugoniots.  The  obvious  conclusion  is  that  quartz  re¬ 
sponds  differently  when  impacted  by  porous  materials  and  cannot  be  used 
indescriminately,  even  if  the  porouc  material  fully  compacts  to  solid 
density  within  tha  rise  time  of  the  shock  profile.  Tha  lower  readings  in 
X-cut  quarts  are  consistent  with  the  piezoelectric  tensor  for  tha  material. 
The  observations  about  forarunner  pressures  in  Section  VI -a  should  not 
bs  affected,  however,  because  adequate  allowance  has  been  made  in  their 
qualitative  nature.  Furthermore,  errors  in  forerunner  magnitude  should 
hav«  tho  asms  (small)  effact  on  the  experimental  points  and  smooth 
curves  of  Figs.  18-20,  eo  that  tha  conclusions  of  this  section  (VI-3)  re¬ 
main  unaffected. 


Shock  Attenuation  Experiments  and  Calculations 


Ona  important  goal  of  die  research  program- was  to  develop  the 
ability  to  predict  the  behavior  of  a  foam  subjected  to  short-duration, 
ht^h -amplitude,  impulsive  loading.  Tha  flyer  plate  experiments  des- 
cribedln  Section  V  were  designed  to  approximate  such  pulses  with 
momentum  densities  of  0.!  to  0,75  x  104  dyne-sec/cm8  (see  Table  IV) 


HUXX  JttUMlKUXV  X3M*  «M  YiVQ 


delivered  within  leas  than  0.  15  /isec.  The  resulting  voltage-time  pro¬ 
files  at  the  quartz-foam  interface  were  recorded. 

If  the  quarts  gage  current  is  proportional  to  the  stress  difference 
between  its  two  faces,  a  reasonably  good  idea  of  the  pressure  profile 
can  be  obtained  even  after  a  shock  has  transited  the  gage,  since  the 
gage  is  backed  up  by  material  of  similar  shock  impedance.  The  shock 
transit  time  in  the  gage  was  1.66  Jlsec,  so  profiles  were  reconstructed 
by  adding  to  each  pressure  reading  the  reading  1. 66  ftsec  ahead  of  it. 
This  procedure  should  result  in  computed  pressures  which  are  slightly 
high,  due  both  to  the  action  of  the  shock  on  the  gage11  and  to  the  fact 
that  the  impedance  of  the  backing  material  was  actually  slightly  lower 
than  that  of  quartz.  On  the  other  hand,  the  observations  reported  in 
Section  VI-3  indicate  that  use  of  the  quartz  gage  on  porous  targets  may 
result  in  erroneously  low  pressure  readings  if  the  standard  calibration 
(Fig.  9)  is  used.  Therefore,  the  reconstructed  experimental  profiles  ' 
(solid  lines  of  Figs.  22  to  24)  should  be  considered  reliable  as  to  timing 
and  general  shape,  but  not  necessarily  as  to  absolute  pressure  level. 

It  should  also  be  noted  that  profiles  at  times  greater  than  about  4  Jisec 
after  impact  (t"0  at  impact)  will  deviate  towards  lower  and  lower 
pressures,  due  to  the  arrival  of  side  rarefactions,  and  should  therefor# 
not  be  compared,  with  the  one -dimensional  computer  calculations. 

To  test  our  ability  to  'predict  the  response  of  the  foams,  calcu¬ 
lations  were  made  for  all  profilas  by  using  ths  artifical  viscosity,  digital 
computer  program  and  equation  of  stata  information  describad  in 
Section  IH.  The  computtt’onal  procedure  sometimes  distorts  the  presw 
sure  at  die  quarts -foam  interface,  in  which  case  the  pressure  in  Zone  4 
of  the  quartz  has  been  used  and  the  times  hevs  been  adjusted  to  those  at 
the  interface  by  using  the  shock  velocity  in  the  quartz  end  noting  that 
each  cell  is  0. 008  cm  vide.  ■*!*•»  resulting  profiles  are  also  plotted  in 
Figs.  22  and  23. 
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FIG.  23  RECONSTRUCTED  AND  COMPUTED  PRESSURE-TIME  PROFILES  AT 
FOAM-QUARTZ  INTERFACE  FOR  SHOCK  ATTENUATE::  STUDIES 
Fhrw  pkrt*  'mputt  accwrs  at  tern  •  0.  P—m  ■*  2.1  “cKuMt"  cImbIrv* 
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FIG.  24  RECONSTRUCTED  AND  COMPUTED  PRESSURE-TIME  PROFILES  AT  FOAM- 
QUARTZ  INTERFACE,  SHOWING  EFFECT  OF  VARIATIONS  IN 
COMPUTATIONAL  PROCEDURE 
F«m  t*  ATJ  jrophlt*,  SW  11,451. 


For  ATJ  g rmphits.  Fig.  22  Indicates  that  agreement  between  cal¬ 
culated  and  experimentally  measured  profile*  i*  very  good.  For  the 
aluminum  foam,  agreement  between  calculated  and  experimental  pro* 
files  is  fairly  good.  Presumably,  agreement  could  be  Improved 
and  certainly  made  more  meaningful  by  a  better  knowledge  of  forerunner 
preseuree  and  any  resulting  modification  in  the  constants  of  Eq.  (17). 

It  is  clear  that  a  porous  eolid  will  evidence  rise  times  that  d«* 
pend  upon  the  scale  of  homogeneity;  i.  e. ,  a  finite  time  is  required  for 
each  particle  of  the  foam  to  come  into  pressure  and  particie-«rvtccity 


equilibrium  with  itc  surrounding*.  On  the  othe-  h end,  the  artificial 
viscosity  used  in  the  computations  acts  like  a  real  viscous  force  against 
which  irreversible  work  must  be  done  in  compressing  the  solid.  By 
appropriately  adjusting  the  value  of  Q,  the  rise  times  of  the  calculated 
pressure  pulses  can  be  changed,  and  there  is  a  corresponding  erosion 
of  the  pulse  because  of  the  additional  energy  dissipation.  With  this  in 
mind,  an  attempt  to  match  the  profiles  in  aluminum  foam  was  made  by 
increasing  the  value  Cona  from  0. 1  to  0.  5  in  Eq.  (11)  and  adjusting  the 
values  of  R  and  K  in  Eq.  (17).  The  resulting  profiles  were- in  better 
agreement  with  experiment  but  the  introduction  of  a  larger  linear  vis¬ 
cosity  as  such  has  not  been  justified  on  a  physical  basis,  further  in¬ 
vestigation  of  the  problem,  especially  from  a  microscopic  point  of 
view,  would  be  of  considerable  interest. 

In  order  to  see  the  effects  of  various  assumptions,  the  experi¬ 
mental  and  computed  profiles  of  Fig.  ZZa  have  been  reproduced  in 
rig.  Z4,  along  with  profiles  calculated  (1)  neglecting  the  styrofoam 
region  behind  die  flyer  plate  and  (Z)  Including  the  styrofoam  region  but 
using  release  paths  for  ATJ  graphite  coincident  with  the  compression 
paths  (NOREL  procedure),  ms  shown  in  Fig.  4  of  Section  HI.  It  may  be 
seen  that  the  styrofoam  had  only  a  small  effect  and  that  die  "spring- 
back"  in  ATJ  graphite  must  occur  at  essentially  sero  pressure,  so  that 
the  "locking"  release  curves  arc  the  appropriate  ones  to  use. 

In  conclusion,  it  might  be  stated  that  while  reasonably  good 
agreement  can  be  achieved,  the  shape  of  the  P-V  Hugonlots  in  the 
partially  compacted  region  has  been  conveniently  chosen  by  Eq.  (17) 
with  appropriate  constants  and  is  sensitive  to  errors  in  assumed 
forerunner  pressures.  It  would  also  be  desirable  to  confirm  ex¬ 
perimentally  the  choice  of  assumed  release  paths.  Another  possible 
source  of  error  lies  in  the  extrapolation  and  normalization  of  data 
for  the  solid  materials .  The  introduction  of  a  linear  viscosity  to 
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aid  in  matching  wave  profile*  can  be  of  value  but  ha*  not  been  rigorously 
justified  on  physical  principle*.  A  possible  refinement  of  the  calculation* 
would  be  to  employ  an  elaetoplastic  model  for  the  Hujoniot*  of  the  com¬ 
pacted  material*. 


SECTION  YU 

SUMMARY  AND  RECOMMENDATIONS  FOR  FURTHER  WORK 

A  portion  of  the  effort  was  directed  to  the  development  of  suitable 
aluminum  foam*  on  which  to  make  shock  measures,  snt*.  Aluminum 
prec ents  a  particularly  formidable  problem  because  of  the  omnipresent, 
nearly  impervious,  oxide  layer.  We  have  succeeded,  however,  in  pro¬ 
ducing  relatively  homogeneous  foams  in  the  density  range  of  40  to  SO 
percent  of  crystal  density.  We  have  also  successfully  varied  particle 
shape,  sixe,  and  size  distribution.  Production  techniques  investigated 
include:  (1)  hot  pressing  of  aluminum  powders;  (2)  cold  pressing  and 
subsequent  sintering  of  aluminum  powders;  and  (3)  repeated  h'>t  pressing 
and  sintering  of  a  mixture  of  aluminum  powders  and  silica  microballoons. 

Hugoniot  measurements  using  newly  developed  techniques  have 
been  made  on  several  graphite  (or  carbon)  and-aluminum  foams.  Hugo- 
niots  of  one  graphite  and  one  aluminum  foam  (both  about  three -fourths 
of  crystal  density)  have  been  traced  in  some  detail  up  to  about  25  kbar. 
Attempts  have  been  made  to  determine  the  pressures  transmitted  to  a 
quartz  (or  solid  aluminum)  structure  for  given  pressures  in  the  foams. 

We  have  found  in  general  that  la  the  pressure  and  porosity  range  studied, 
the  "compacted"  volumes  for  pressures  above  a  few  kilobars  are  essen¬ 
tially  those  of  the  solid  materials  at  the  same  pressure . 

In  the  case  of  aluminum,  (ue  foam  appears  to  "lock"  to  approxi¬ 
mately  solid  aluminum  density  and  lose  Its  fosm  characteristics  (i.  *. , 
the  density  of  recovered  shocked  fosm  corresponds  approximately  to 
that  of  solid  aluminum).  On  the  other  hand,  recovered  shocked  speci¬ 
mens  of  ATJ  graphite  exhibit  densities  very  close  to  their  initial  den¬ 
sities  (!•••»  they  are  recovered  as  foams).  The  equation  of  state  in¬ 
formation  outlined  above,  is  of  rmjor  Importance  for  any  shock  attenu¬ 
ation  calculations . 
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,  Parameter  variation  studies  included  various  particle  (or  pore) 

•ices,  shapes,  and  distributions.  In  general  the  preliminary  results 
indicated  that  these  factors  are  probably  of  less  importance  than  the 
strength  of  the  solid  matrix.  Where  the  matrix  is  strongest  the  fore- 
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runner  speed  and  amplitude  may  be  expected  to  be  the  greatest.  There 
is  a  tendency  for  these  desirable  characteristics  to  be  associated  witk 
foams  of  higher  density  if  the  other  factors  are  kept  constant.  With  the 
present  state  of  material  development  of  foams  of  the  type  studied, 
there  is  too  much  scatter  in  forerunner  characteristics  to  enable  more 
specific  correlations. 

The  artificial  viscosity  computer  code  has  been  successfully 
adapted  to  calculation  of  shock  attenuation  in  porous  solids.  Within 
the  idealisations  of  the  models  employed,  calculated  transit  times  and 
■hock  profiles  are  in  reasonably  good  agreement  with  the  experimentally 
measured  quantities,  particularly  for  the  case  of  ATJ  graphite.  However, 
uncertainties  in  forefunner  pressure  and  in  the  Hugoniot  curve  in  the 
partially  compacted  region  limit  thir  reault  to  a  demonstration  of  capa¬ 
bility  rather  than  providing  a  definite  numerical  result. 

During  tha  course  of  the  work  several  topics  have  suggested 
themselves  ae  Important  areas  for  further  research: 

1.  Precise  measurement  of  forerunner  pressures  (e.g. ,  by 
using  an  aluminum  witneas  plate  between  the  foam  and  a  tarts  gage) 
and  of  pressure-time  profiles  for  attenuation  shots  should  be  mrde  so 
that  thr  actual  attenuation  profiles  can  be  established  and  checked 
against  the  artificial  viscosity  calculation* . 

2.  Research  on  new  materials  such  as  heavy  metal  foams  and 
ceramic  foams  should  be  initiated. 

3.  The  study  of  pertinent  structural  variables  controlling  the 
forerunner  speed  and  amplitude  should  be  extended  and  expanded. 
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Correlation*  with  various  static  and  quasi-static  measurements  should 
continue  to  be  investigated, 

4.  Optimum  distention  ratios  for  a  given  material  should  be 
established.  This  optimum  ratio  will  depend  upon  the  type  of  foam, 
but  for  a  given  foam  it  will  involve  a  compromise  between  the  large 
forerunner  amplitude  (and  speed)  associated  with  low  distentions  and 
the  rapid  "attenuation  from  the  rear"  associated  with  higher  distensions. 

5.  Hugoniot  measurements  should  be  extended  to  wider  porosity 
ranges,  and  an  attempt  should  be  made  to  obtain  reliable  data  at  pres¬ 
sures  just  greater  than  that  of  the  forerunner  wave.  A  knowledge  of  the 
actual  pressures  required  for  complete  compaction  is  very  important 
for  the  calculations  of  shock  attenuation.  At  present  we  can  assign 
upper  and  lower  bounds  to  these  values  for  graphite  and  aluminum  but 
we  do  not  have  a  direct  measure  of  the  P-V  curve  between  elastic  and 
compacted  segments. 

6.  Hugoniot*  for  aluminum  foams  with  different  particle  sixes 
and  shapes  have  been  found  to  essentially  coincide  at  about  4  kbar,  hut 
the  work  should  be  extended  to  other  material*  and  to  higher  pressures 
where  sny  differences  rrAght  be  more  easily  observed. 

7.  The  release  curves,  especially  in  partially  compacted  re¬ 
gions,  are  of  major  importance.  It  is  also  clear  that  the  release 
Curves  for  compacted  graphite  are  not  fcimplc  locked  curves,  but  the 
way  in  which  they  deviate  from  the  locked  model  has  only  been  inferred, 
and  has  not  been  experimentally  measured.  Attempts  should  be  made 
to  measure  release  curves  directly. 

S.  A  basic  inquiry  as  to  why  some  foams  (like  aluminum)  essen¬ 
tially  lock,  and  others  (like  graphite)  do  not.  would  be  of  considerable 
interest.  It  would  further  be  of  interest  to  discover  whether  a  shocked 
foam  which  dues  not  lock  may  act  at  an  effective  countermeasure 
malarial  in  case  of  subsequent  "hits.  "  Investigation  of  release  curves 
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may  shed  some  light,  but  experiments  designed  so  that  specimens  with 
a  well  known  history  can  be  recovered  intact  for  detailed  study  and  per¬ 
haps  subsequent  shocking  would  be  valuable. 

9.  Accurate  Hugoniots  for  the  solid  materials  should  be 
measured  at  low  pressures  in  order  to  eliminate  the  uncertainties 
associated  with  extrapolation  and  normalization  of  data.  The  elastic- 
pla  Stic  model  should  then  be  used  in  the  computer  calculations. 

10.  Investigation  of  dissipative  processes  in  foams  should  be 
conducted  from  a  microscopic  viewpoint  to  lead  to  a  better  understand¬ 
ing  of  wave  profiles. 


APPENDIX  I 


DESCRIPTION  OF  CALCULATION  ROUTINE 
USING  THE  ARTIFICIAL  VISCOSITY  CODE 

A.  Essential  Step* 

The  essential  steps  in  the  calculation  are  summarized  below: 

1.  Read  in:  Material  constants,  number  of  regions,  dimension 
and  kind  of  material  in  each  region,  number  of  cells  in  each 
region,  geometry  (plane,-  cylindrical,  or  spherical),  Q -con¬ 
stants,  control  instructions  for  printing  output,  flyer  plate 
velocity. 

2.  Determine  cell  size  in  each  region. 

3.  Set  initial  conditions  for  each  cell:  mass  in  each' cell, 

Euler  coordinate  of  each  cell.  -  - . 

4.  Begin  computation. 

5.  Calculate  new  particle  velocity  for  each  cell  from  momentum 
equation.  * 

6.  Calculate  new  Euler  coordinate  for  each  cell. 

7.  Calculate  new  specific  volume  from. continuity  equation. 

8.  Calculate  new  Q  for  each  cell.  ^ 

.  Calculate  new  F  for  each,  cell.  .  ^  ,V~' 

10.  Calculate  At  for  next  time  etep.  4 

11.  Print  result*. 

12.  Return  to  step  5  and  repeat.  Continue  cycling  until  proUem. 

5  .  is  complete;  x-  ■  ;  .  , 

Both  the  initializing  (Step  3)  and  the  computation  of  P  (Step  9}  are 
accomplished  in  procedures.  The  initialising  is  straightforward;'  the  ’  , 

computation  of  P  requires  iteration  if  the  constitutive  relation  is  energy- 
dapA-rient,  If  it  is  not  energy -dependent,  it  is  computed  directly  from 
the  specific  S'olumt  previously  determined.  The  entire  program  ae  used 
in  a  typical  run  ie  reproduced  at  the  end  «f  this  appendix. 
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B.  Zoning 

The  experimental  configuration  i*  described  in  Section  V.  Brief¬ 
ly,  it  i*  a  one  -dimenaional  experiment  which  consists  of  a  flyer  plate 
assembly  (a  6061-T6  aluminum  flyer  plate  backed  by  a  2  lb/ft3  styro¬ 
foam  support)  impacting  the  target  (foam)  which  is  bounded  on  the 
opposite  side  by  quarts. 

Computation  is  started  at  the  instant  the  flyer  plate  strikes  the 
target.  Initial  values  are:  pressure  PJ0  everywhere;  particle 
velocity  u  *  u^*  constant  in  flyer  plate  assembly;  u  -  0  elsewhere; 
specific  volume  has  its  zero  pressure  value,  V  *  V0»  everywhere. 
Artificial  viscosity  Q  is  zero  everywhere.  A  value  of  the  time  incre¬ 
ment  At  (DELT)*  must  be -specified  in  the  input  for  the  first  cycle. 

Care  must  be  taken  that  this  not  be  too  large.  If  it  is.  a  pressure  spike 
is  produced  at  the  flyer  plate/ specimen  interface,  which  causey  a  wildly 
fluctuating  output. 

The  number  of  zonae  in  «ach  region  can  he  specified  somewhat 
at  will;  however,  fluctuations  in  lhe  output  are  reduced  by  choosing 
zone  size  for  each  region  in  such  a  way  that  the  minimum  At  (TLIMA) 
for  each  region  has  about  the  same  value.  Zone  size  ie  calculated  as 
initial  region  thickness  divided  by  number  of  zones  in  region  (under 
SAMES);  here,  too,  reference  is  made  to  procedures  ZSOLID  and 
ZWATERr  in  which  initial  values  are  set.  . ; 

Geometry  is  msde  planar  by  setting  ALP  *  1  under  LXSTMTLS 

A  logical  variable  is  set  PUP  [j]«  TRUE;  its  use  is  described 
In  Section  *>  of  thie  appendix.  -  '  1. 

C.  Computing  Steps  - 

The  computation  for' a  given  cycle  starts  at  LEFT:  PPEAK*"  0, 
etc.  At  this  point  all  parameters  for  the  previous  cycle  (time  »  t) 
have  been  determined,  t  bat  been  advanced  to  *  +  At,  and  cycle  to  cycle 
+  i.  •  • 


*  Labels  amd  names  of  variables  and  procedures  have  reference  to  the 
ALGOL  program  used  for  the  computations. 
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Computation*  are  made  in  the  following  order: 

Utl]  i*  a  new  value  for  particle  velocity  in  the  first  zone, 

J  «  1. 

XA  is  temporary  storage  for  the  value  of  the_  Euler  coordinate 
of  the  Jth  son*. 

Ufjl  3  is  the  new  value  of  U  for  the  (J  +  1)  zone. 

Xfj]  is  the  Euler  coordinate. 

VN  is  the  new  value  of  specific  volume  for  zone  J. 

QA  is  the  new  value  for  stress  from  artificial  viscosity. 

New  values  of  P[j3,  ECj3  are  calculated  in  the  procedures 
styrofoam,  ELASTOPLASTIC,  WATER,  or  PHASETRANSITION;  the 
procedure  entered  at  this  point  depends  upon  the  value  of  the  parameter 
BURNTsI,  which  has  previously  been  prescribed  for  each  of  the  four 
regions:  S  *  1  (styrofoam  support),  S  •  Z  (aluminum  flyer  plate), 

S  *  3  (target),  S  *  4  (quartz);  i.  ♦*. ,  • 

3URN[l3  »  5 
BURN [2]  «  0 
BURN[3]  «  4 
BURNt4]  »  3 

Following  the  calculation  of  P[j],  the  variable  V[j3  is  replaced 
by  the  number  temporarily  stored  in  VN,  and  J  is  tested  to  see  whether 
it  is  beyond  the  shock  front.  If  it  is  not  (FALSE),  another  J-cycle  is 
made  in  response  to  the  instruction  "IF  J  *JSTAR+1  THEN  GO  TO 
GETUV. "  If  It  is  (TRUE),  a  new  At  is  determined,-  the  information 
stored  for  the  past  cycle  is  printed  (if  desired),  and  a  new  time  cycle 
is  started  by  the  last  instruction  of  the  program,  "GO  TO  LEFT. " 

D.  Comments  on  PHASETRANSITION  Procedure 

The  implementation  of  the  prescription  for  calculating  P  and 
V  (see  Section  Hl-D)  is  simple: 

Af  «  V  (J,t*At)  -  V(J,t) 

can  be  stored,  and  as  long  as  it  i*  'negative,  compression  of  that  zone 
is  occurring,-  When  it  turns  positive,  expansion  has  started, 
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Although  we  have  not  found  it  to  be  the  cate  to  far,  it  is  possible 
that  fluctuation*  inherent  in  the  numerical  integration  may  eau«e  diffi¬ 
culties  with  thi*  simple  criterion.  To  protect  against  that  eventuality 
if  it  should  arise  in  the  future,  a  routine  has  been  included  in  the  pro¬ 
gram  to  smooth  out  the  effects  of  fluctuations  in  the  sign  of  AV.  The 
disadvantage  of  its  use  is  the  resulting  delay  in  changes  between  com- 
pression*and  rarefaction,  so  it  is  recommended  that  the  value  of  n 
(see  below)  be  kept  as  small  as  possible.  We  have  used  m=l  (which 
corresponds  to  the  simple  criterion  above)  in  all  cases. 


Routine  to  Check  for  Compression  or  Rarefaction.  Running, 
average  values  of  V  and  of  Av  could  be  carried  along  and  these  should 
have  sufficient  stability  to  serve  as  controls.  The  normal  definition  of 
a  running  average  over  an  interval  (t-T)  to  t  can  he  written 

VR(t,T)«  ^ 


vNtjj.m) 


N 

kl.  V 

n  *  N  -  m 


v(tA> 


VR(tN,m)  «  VR(tN-l)  +  V<  W*  ]  (26) 


If  this  were  carried  in  the  present  computation,  it  would  require  that 
values  of  PCl]  for  the  previous  m  cycles  be  stored.  In  order  to  con- 
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If  we  denote. 
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V(t^)  by  Vjj7~etc. ,  £q*>  (26)  and  (27)  can  be  written 

(28) 


7*  -  VR 


VN-VN'm 


and 


N  1  N-l 


'N  'N-l 


m 


'N 


rN-l 


m 


(29) 


Then  if  V  passe*  through  a  minimum,  VR  passes  through  a  minimum 
somewhat  later.  But,  unless  changes  have  been  very  rapid, 

Vl<VN-  m;  so  7  passes  through  its  minimum  after  V  and  before  7  R. 
Consequently,  it  should  serve  as  an  adequate  monitor  for  changes  be¬ 
tween  compression  and  rarefaction  if  m  is  chosen  properly. 

The  monitoring  function  of  7  is  fulfilled  by  using  V  to  determine 

a  reference  "turnaround"  volume,  V  ,  such  that  whenever  V  >  V  ,  the 

c  c 

state  point  is  assumed  to  lie  on  a  reversible  curve.  CF  of  Fig.  3;  and 
the  function  describing  these  carves  is  f^f(V),  given  by  Eos.  (18)  and 
(20).  The  initial  elastic  curve,  AB,  of  Fig.  3,  is  thus  included  in  f^ 
(V).  If  V  it  ever  *  it  is  assumed  that  the  material  is  locked  into 
its  solid  state,  and  all  subsequent  compression  and  expansion  will  take 
place  along  the  curve  f  j(V)  given  by  Eq.  (16).  Thus,  for  each  cycle 
the  function  e£  the  machihe  is  to;. 

1)  Calculate  new  7  and  A 7*  ; 

t  * 

2)  V  new  7<VC,  then  calculate  P  *  fj(V)  or  f2(V)  or  fj(V) 
and  return  to  main  program 

3)  tf  new  V>VC,  then  set  VC.  PC  and  VCO;  calculate  P  * 
fj^tV);  and  return  te  main  program. 

The  information  stored  in  the  memory  at  the  time  of  catering 
the. procedure,  (with. definition  of  symbols)  is: 

PUPtj]  «  TRUE  or  FAgJSE.  If  TRUE,  P  »  fj(V)  cr  f2(V)  or 
f,(Y).  If  *AL5*,  fc«ftot(V).  ~ 


n. 


£  * 


.  A 


-^•T- 


T  \ 


INT[j]  3  TRUE  or  FALSE.  If  TRUE,  VC,  PC  and  VCO  have 
not  been  cet.  If  FALSE,  they  have  been  aet. 

PTj]  *  P(J,N-1}  *  previous  value  of  P 
VBArCj}  *  ^(J,N-2) 

PCLJ ]  3  PC(J) 

VClj]  =  VC(J) 
vcofj]  *  Vco(J) 

MM  *  m  *  number  of  valuea  of  V  in  P 

VN  3  V(J,  N),  the  current  value  of  V  to  be  used  in  computing  the 
•new  preaaure. 

Parameters  to  be  determined  in  the  procedure  are: 

VBARA  3  F {J, N-l)  and  VCtlX  PCCj],  VCO[j]if  P  turns 
around 

PW  3  new  value  of  P  3  P(J,N) 

.  pupCj!  ...•••  ' 

INTM  • 

E.  Comment  a  on  NOKEL  Procedure 


The  NOREL  procedure  if  very  simple.  A  flow  chart  is  shown 
in  Fife.  25,  since  the  procedure  does  not  appear  in  the  program  listing 
in  Section  T  of  Oils  appendix. 
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BE6IN 

FILE  KIN  (3.10).  POUT  *M3»15)» 

INTEtoEK  NUN)  LABEL  FINISH! 

R£AD(R1N»/»NUH) I 
NUM  *  NUN  ♦  101 
BEGIN 

INTEGER  COUNT . COUNTS. Clr  CLE. C YCLES » HALT » J.J1 'UC.JCC.JCRIT# UMAX .RESTART » 
JPMAX.USTAR.UT . NT. NTT. OPTION. S» SO » SX  »EF» 

SP 1 1 »  JPB » JPE  #  ALP  #  RR I 

REAL  CONA  #  CO  *  CttStt  *  CQSQ4 . DELT  •  OELTIr  OELU  .  DELX  »  DfNU.DL. 

OLi . OH . DPOMU . DTMX . OTN * OTNH1 . QA » JMAXF » 

LcFTP .LINEAR. LX »M# Ml »MUN.P£X»PF.PH#PLEFT  »PMAX*PMIN»PPEAK» 
PSCALE.TAU. TIMES# TLIMB»TOUXT#UO» 

2>Pm.EW.F«»EA.EU. VJ.BB'TEST.Pll. Vll.KK. 

CsPS.CSHE# VN.XO. YE1# YF»  YPM.XA.XSPA.VR  I 
INTEGER  ARKAY  CELLC0:253#eURN.Ht0J9H 

REAL  ARRAY  A.8.C.DX.6.KUA.MY.PA.PKA.RH0. VO. Y.OV.GAM. VCJ.LE0153. 

AAt  o:b3.Ttf(8:253»£»F.MUt..MUF.P.PEf O.TMP. EI‘?I .X* 
XALPS.ENT.CSP.TB.TLIMA.U.YP.  YC0SNUN3I 
BOOLEAN  BOOL I 

COMMENT  PHASETKANSITION  VARIABLES I INTEGER  MM » REAL  Vl.Pl.V2.P2. V0.V6.K.RI 
REAL  .yBAHA.DVBARAI 

REAL  AHRAY  VBAK.0V8AR.VC. PC. VCQCOSNUHJI 
BOOLEAN  AKRAY  PUP»1NTC0:NUM3I 

FORMAT  FHTD(X10»«TIM£  . CIA »B*X4. "CELT  s".EHu®.X«»»"DTN  :*.E1«.8#X*# 

"CYCLE  =".I5#X4»"UCRIT  =".I5./>»  • 

0TAG<X2."J*»Xa."U"»xn."V»»XXX."P"»Xn»*E"»XU."Q"»Xll» 
"X"»X10."NUe"yX9."NUF"«X3»"TUMA"»//>»  . 

F3(X3.0F12.6.E13.S>  •  -- 

Fl(X6."ALP".X10»"DeLT"(Xll.«OTMX".XXl»"CONA"»XX3."CG"./.I8.X7. 
AE1S «*.//)» 

F*( "CYCLES  COUNTS  NTT".XG»"TtrC  3. •"</» 3I6«X2.10F10«3./»10(X20* 
X0F10.3*/)*/)*  . 

FBUHN  (110. "REGIONS*. 915#/)* 

?•  (X8."A".Xll»"B".XH»"C".Xn,"G".XX0."MUA".XX0»*PA**X10."HY". 

AlO.'*RHO"»Xi  »*H",»Xll»"L"). 

FI  (8F12.4»ia»X4*  FA2.4.//1. 

HPmX4»"N*»X5."a*»Xt."Vi«»XU»«PI«.m."V2»*Xil»"P2"»XiX»,,V0"»Xll 
»"VCf‘.Xll»"K“.Xl0»"H"»X5»"L"»/»2l6»7F13.6.l6.Fll»6»//)» 
Hti<X7»"¥r0*»X9»*H*»XS»"L,'»/»Fl3oG»lB»Fll.S»//l. 
FPUt//.X10.*TAU*»Xl5»"LEFTP*»X16»"UCl  J".X9»<'0PTI0N,,»/» 
3ESQ.6.I3.//)* 

FSUM-TIME  s*»ll%.a»X4."UFS  s»»EX4.e*XG»nXFS  =*,.E1*.«»X%» 

"CYCLE  S"»IAtX3."JPMAX  =".IG.X3.*JCRIT  s"»IWI» 
F(*AGE(U.BFXS.»)» 

LIST  INTERVAL! TlMES.OELT.OTN. CYCLE. JCRITJ. 

EMU«I.UIj3.X(J}Ii 

MTLS  (ALP.OELT.OTMXyCONA.CO)  r 

CONTROL (CYCLES»COUNTS*MTT» FOR  NT*X  STEP  X  UNTIL  NTT  00  T*CNT3>» 
STREAMLINE! JCC. FOR  JC»1  STEP  X  UNTIL  JCC  00  CELUJC3). 
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RtGIONS  (SI • FOR  S  *  1  STEP  1  UNTIL  SI  DO  BURNCS3) 

£QST<ACS3»6CS3«CtS3.GtSJ*MUACS3#PACS3»MYCS3*RHOCS3#HCSj»LtSa>* 

PHAS£(MM»R.Vl»Pl»V2#PS!*V0»VS#KtHCS3*LCS3)# 

QUAKTZ(V0CS3»HCS3»  LCS  j)  * 

PULSE ( TAU*LEFTPrUC 1 3*  OPTION) * 

ARRAYS*  J*UCJ1*  VC  J3*PC  J3#Et  J3*QC  J3.XC  J3*KUECd3*MUFC  J3#TLIMA£  J3)» 
FACE(TXM£S#Ut  J3»XC  J3»CYCLE»  JPMAX*  dCRXT)  t 
FORGAGEt JiUCU+1  3->VC«.ID*PC  J3*QC  J3»XCU*13)  I 
LABEL  KEaOING»L£FT*SETJ.ULEFT«6ETUV»GETX.FLAP.NEWT) 


PROCEDURE  /WATER  I  BEGIN 
CSPS  ♦  0.3  I 
BB  ♦  0.11 

FOR  J  *  (HC5-13  ♦  1)  STEP  1  UNTIL  HCS3  DC  BEGIN 
VI J3  *■  VOCS3  I 
Xld*13  •  XCd3  ♦  OXCS3  ) 

KALPSCJ3  ♦CXCJna  *  ALP  -  XCd3  *  ALP)  /  VCJJ  I 
ut j*i j  ♦  eCJ3  *  o.«  i 
PCJ3  *  0.01 
TLIHAL J3  *  DECT  ? 

CSPCJ3  •  CSPS  * 

VBARC J3:~  VO) 

PUPCJ3  ss  INTCJ3  :s  TRUE! 

.  tNyl  ENO  /WATER  I  • 

PROCEDURE  WATER) 

COMMENT  QUARTZ  REPLACES  WATER)  - 
BEGIN 

PtJ3  •  .OG9:>kU.O-2.6S*VN)I  - 

CSPCJ3  «•  .5720) 

8CU3  •  «A)  :  . 

•ENO  WATER)-  * . rv  ‘  \  ’ 

PROCEDURE  STYROFOAM)  , 

BEGIN 

IF  VN>31»8S  THEN  BEGIN  VN«31«*5)  PW  *  0.0)  END)  i 
1F  31.252YN  AND  VN>. 942*1  THEN  PW*<3l.25-VN)/<31.25-.94221)XO-3 
*(1.0  -.lx(VN-.94221)/(31.25-. 942*1))) 

IF  VNS.94221  THEN  PW  »(. 390470  ♦.*0fll47x(.95*4/VN  -1))*(.9324/WN  -1)» 
KdJ*Pw)  .• 

«CJ>GArrt'V# .  .  -  •  ■  ■ 

CSPtd)  .4)  v-'-' 

ENO  STYROFOAM)  ; 

PROCEDURE  /SOLID  )  BEGIN 

VOtSJ  *  1.0  /  RHXS3  ) 

CSPS  ♦  SORT  (0.0  ♦  4.0  x  0CS3  /  3.0  *  V0CS3  x  ACS3)  t 
Fat  d  *  (HtS-l)  ♦  U  STEP  1  UNTIL  HCS3  00  BEGIN 
VtdJ  ♦  VOCSI  » 
xtduj  •  xtdi  ♦  ottsa  i  ■ 

XALPSCdl  MXtdMJ  •  ALP  -  Xtd3  •  ALP)  /  VCdl  »••  ; 

.  Utd+13  *  PCd3  ♦  Xd3  Ml 
MjFtdS  *  -0.0001  ) 

MUEXd)  *  MUACSa  > 
pad)  *•  PHACS)  I 
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YPCJ3  ••  YCSJ  * 

TLIKALJJ  DELT  l  ■ 

CSPtJJ  ♦  CSPS  I  •  - 

EL  J  J  *■  0.0  » 

END  EMU  Z50LID  * 

S 

PROCtOURc  PHASfcTRANSXTION* 

BESIN 

FORMAT  F7(9(F11.6«X2>>* 

LABEL  RESET I 

VBARA: = ( ( MM-l ) xVBARC J 3* VN> /MM) 

UVBAHA:=VBARA  -  V8ARCJ3I 

IF  NOT  PUPCJ3  THEN  IF  VNSVCEJ3  OR  VCtJV.y  THEN  PUPEJ3  *  INTXJ3  *  TRUE* 
if  PUPLJJ  THEN  If  VBARA  >  V6ARIJ3  THEN  P,.$P£J3  »  FALSE* 

IF  PUPLJJ  THEN  BESIN 

IF  VI  SVN  THEN  BESIN 

Pur’-PIMI.O+BBxXVN-Vl  >/<Vl  -VO) )*( VN-V0)/(V1  -V0)» 

GO  TO  RESET*  END* 

IF  VNl  V2  THEN  BESIN 

Pi  *  PI*  K«((V1/VN)*R  -1»/R* 

GO  TO  RESET*  END* 

END  I 

IF  VfKVZ  THEN  BESIN  PW*(.192?2  ♦2.85475«XVS/VN-ll )*<VS/VN>1) » 

COMMENT  THIS  EQUATION  IS  FOR  GRAPHITE*  FOR  ALUMINUM  SUBSTITUTE 
P*  *  (.72922  ♦  1.71742* (VS/VN  -1)}k{VS/VN  -1>» 

60  TO  RESET *  END* 

IF  (0VBARC«j3<0  ANO  2NTLJ3)  OR  VJSV2  TKiiN  BESIN  PCC J3:=PC J3)  VCC<J3:=VCJ3* 
IF  VJCV2  THEN  VttJl  ♦  V2* 

VC0CJ3  *  VO  ♦lVG-V0»*tVCCJ3-Vl»/CV2-Vl>*  IMTCJ3  S=  FALSE* 

ENO* 

IF  VCIU3-VC0C J3  THEN 

BEGIN  vAlT£tP0UT»F7» J»VCC J3»VC0t J3* VN) ISO  TO  RESETtENOl 
IF  VIOVCOCJ3  THEN  VN  ♦  VC9CJ3* 

PB»PCEJ3*U^(VN-VC£J3**ClF.l*(Vf»-VCOCJ3)/(VCEJ3-VCI)CJ3*»/<VCtJ3-VCOC4 

i»M 

IF  P«<0  THEN  BESIN  P*+*.9 *  VN  ♦  VC0CJ3*  END*.  - 
RESET!  - 

IF  ABS<VN-VJ)>0-3  THEN  BESIN  -  • 

TEST  ►  ~(VN*2l*«P*f-PtJ31/<VN-VU»f  - 
IF  TEST  2  S  THEN  CSPCU3  *  S«RT<TEST>*  ENO* 

9BARCU3  »  VBARA*  KJ2  -  PS*  S£U3»0A*  0VBARCJ3  !=  C VBARA* 

ENO  PHASETRANSITION* 

PROCEDURE  ELASTOPLASTIC* 

BESIN 

LABEL  GETM.CALCM. RETURN I 
StJ3  *  _SAt 
MUN*IV0CS3/VN'**»1»8> 

OPUNU*-MUN*  (  2  .  SKBC  S  a«3.  SxCC  S  3XMUN  t  ♦«[  S  3 1 
CSPC J3*S0RT ( U.8+4. I*SC S 3/3.0 ) xVOC S 3*0P0NU) * 

Pm  MUN*  (  HUN*  l  CC  S  JMMUN*8C  S  3  *  ARC  S  3  »  * 

IF  IWUNJHUEt J3)  AH)  CMUN2MUFC J3I  TH5N  6CSIN 
Pi  J  3*  1 2 . 0*YPC  J  3-4.  t*SC  S  3*PEt  J  3+PHM  4.  0*8CS  3+3. 0»  )/3.S* 

BO  TO  RETURN  ENO) 

TPmNYt  S  3*  tPH-PHAt  3  31  nt  S3* 
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iF  MUN<MUFLJ3  THEN  BEGIN 

MUFC J3*M1*MUNI  M«-O.Oa+Ml*  CL1«— TPH/6CS3* 

6ETM:  PEI»Mx(Mx(CtS3xM«-B£33J*AES3>* 
r£l*KrtS3x(PEl-PHACS3)*US3l 
S3U— t  ( YE14  YPH  > /BL  S  3)  *0 . 5-P*+H»El  I 
OMM (M1-M1/(DLI-DU IxQL* 

Ml  'Ml  DL1-OLI  M»-M-DMI 

IF  ABS ( DM/Mi I >1 • OB-5  THEM  BO  TO  6ETM* 

HUECuJ*MI  PC J3* (-2.0XYPH/3. OI+PHI 
PEL JJ»PE1*  YP£ J3*YE1I  60  TC  RETURN  ENOI 
COMMENT  MUM  IS  GREATER  THAN  MUEIJ3.  SET  VALUE  OF  MUFCU3* 
MUtlJ3*-Ml*MUNI  PEE.  J3*PHJ  YP£ J3*YPHI 
Plj3*-l2.jDxT*»H/3.0M-PH*  M^-O.Pl+Ml*  DL1'>YP,1/BCS3» 

CALCMJ  PF»-MxlMx<CCS7!*M+8tS3)‘>AC£3** 

YF  -M  Yt  S  Jx  (  PF-f  'MAE  S  J  >  ♦  Yt  S3* 

DL-0 . 5x  ( TPH+Yf  I  /6C  S3+PF-PM* 

OM- ( ( M 1-M I / ( 0L1-0U ) xOt* 

Ml-M*  Oul*OL*  M*M-DMI 

IF  ABS <UMFMi»l. OB-5  THEN  *0  TO  CALCMI 

MUFl  J3*MI 

RETURNS  END  ELASTOPLASTJC* 


RESTART  *0*  BOOL  *  FALSE* 

RtA0(RIN./.MTLS)CFINISH3* 

HEAD  IN(*: 

IF  BOOL  THEN  BE6IN  ALP:=AM.i3;OELTJ=AAC23»DT*«:=AA£53lCONATSAAtO» 
CQ:=AAt53*ENOl 

WRITE  (POUUPA6E3)*  WRITE (P0UT*F1»KTLS  J? 

READ C RIN» /» CONTROL* * 

WRITE(POUT«FK.CONTR0U*  •:  .’>•••. 

READIRINt/. STREAMLINE) I 

COMMENT  NT  MUST  BE  SET  TO  l  SO  RESTART  WILL  WORK  »  ■ 

NT  *  l  * 

COMMENT  IF  RESTART  s  1  BO  TO  LEFT  FOR  CONTINUATION  OF  RUN  * 

READ  IRlN*/»RE«IONS*  * 

WRITE  lPOUTCOBL3*FBURN»RE*IONSJ  1 
XCX3  *  0  * 

•  «m  *  o  *  *  - 

BEGIN  LABEL  BOTTCM*  SAMES#  LAST#  NEXT*;..  \  -  •=  *=.  *• 

FOR  S  *  1  STEF  1  UNTIL  SI  00  8EBIN 

sames: 

IP  WURNLSJ  s  0  THEN  BE6IM  ■  > 

R£AD1RIN»/»E*ST> J  ’  r  s  ’’  .  . 

OJttSI  *  US3  /  HCS3 
KS3  *  KS3  ♦  KS-13* 

COMMENT  PA  IS  STRESS  AT  YIELD.  PHA  IS  STRESS  ON  HYDROSTATIC  CURVE  AT 
SAME  STRAIN  AS  FOR  PA* 

YIS>X.««etS3xPACSVU.»+*.QM6CSV5.0»*  HMtSW.SXiftSIrtml 
WMITEIPOUT  »FAI * 

WRITCIPOUT »FB  *  EtST** 

ESOLtO  *  66  T9  BOTTOM  *  END  A 
3 *  BURNCS)  -  *  THEN  BEBIN 


»i 


HEAD  IR1N./.  8UARTZJ*  • 

0XIS3  *  Lt S3  /  HT.S3  * 

HLS3  *  MCS3  +  HCS-131 
tf R 11 E t POUT  » MO  *  QUARTZ* t 
ZWATER  I  GO  TO  BOTTOM  I  END  I 
IF  OURNtSJ  =  M  THEN  BEGIN 
READ  (K£N./«  PHASE! 1 
VOtS3*sVOJ 

DXCS3  *  LCS3  /  HtSI  * 

HI  S3  *  HCS3  ♦  M£ S-l  J» 

KKITEIPOUT*  «PH»  PHASE! I 
£*AT£R  »  GO  TO  BOTTOM  *  £NO  ! ■ 

IF  dURN£S3  =  5  THEN  ‘BEGIN 
READ(RIN'/rLtSJ»HCS3)« 

OXIS3  *  LCS3/KCSI* 

Ht  53  *  HCS3  *  HI  S-l  31 
wOtS3  *•  51*251 
FTATER3  GO  TO  BOTTOM*  Ett»l 
BOTTOM  :  END* 

VtJ3  *  PIJ3  *  01 J3  *  NUEIJ3  ♦  NUFI J3  *  O  t 
READ  (RXSUAPULSE*  i 
WRI TE tPOUTtTPO  »PUUSE1 * 

A  READ IRIN*/-»F0R  I'-X  STEP  1  UNTIL  5  00  AAE13)tLAST3* 

*  80  TO  NEXT* 

LAST*.  BOOLtsTRUC* 

NEXT* 

CLOSE  *R2N»RELEASEl  * 

FOR  S  *  1  STEP  l  UNTIL  SI  DO  IF  BURNT  S  3=8  THEN  EF*HCS3l 
IF  1X13  X  0  THEN  BEGIN  J5TAR  *  HC13  ♦*  I 
FOR  <J  *  1  STEP  1  UNTIL  EF  00  BESIK 
Ut -HI  3  *  UU3  *  END  *  UCJ3  ♦  I.S  «  UCJ3  I  END  : 

ELSE  JWXN 

IF  ALP-*  1  THEN.  ' 

UC13  0X113  i  TSAHCT3  ♦  3..0*  *  ,JSTAR  -*  S  I  END* 
CURRENT  IF  INTERFACE  VELOCITY  DIFFERS  FROM  PLATE  JfELOCXTT  SET  UCJJ  TO 
PSOPER  VALUE  HERE* 

NRITE*P8UT#BTAS)« 

FOR  J  *  1  STEP  1  UNTIL  XHC9»»  DO 
NRITEIPOUT  •■FStARRATSI  * 

FOR  S*1  STEP  1  UNTIL  SI  DO 

FOR  J  *  tHW3~l*  STEP  4  UNTIL  ItitS**)  00 

0RITEtP0UT»F3*A9^AT«»  „  -  - 

WRITfctPCUTtPMtCIt  * 

C8S0*C0«t*  C&r<!"v**.8«C0S»l  LINEAR*!  .Q*CWtA*CONA*  . 
MNCS*CTCL£*«»CRIT*  3» 

OWHOELT*  netn^WlLfMOELT*  TQLttT*A6SCTMCMYT3M  . 

CNO*  •  - 

LEFT*  PPSA**»<  tUNPI.il  TIMES*  TZHE£*D~LT*  «TCLE*CTOE*l* 

setj*  <a*s*i*  ji*«hi*  ^r*uin«  ploft^oi 

CLEFT  S  IF  ALP  *  1  THEN  OEM* 

IX13*»*ttPtI>HEl3-rtJGFf>/IX£t3»Xtl3J  J*"TIl3wOELT2*(jCi3*  END  * 

KA  ♦  Utl3  *  BELT  ♦  *tl3  * 


GETUV 


IF  J  >  HCS3  THEN  S  *  S  ♦  1  I 
OEHU*- 1 XC  JT  J-XC  JU )  /VC  J1 3*  <  XC  J1 3-XC  J  3  J  /VC  J  31 
Ut  U1  JMDELTIx(  PC  J  3-PC  J1 3+QC  J3-OC  J1  J)  >  /OSNUHIC  J1 31  ' 


tfETX  i 

Xlt<3  *  XA  I 

XA  *  CELT  x  UCJ13  *  XtJlJ  J 

IF  J  s  HCS13  THEN  XC  J13  *  XA  I 

IF  AaS(UlJ13)<5.08-5  THEN  UT J13«t)l 

VN  *  (XA  *  ALP  -  XCJ3  *  At?:-  /  XAUP5CJ3  1 

DELU  *  UCJ13  -  UCJ3  I 

OELX  XA  -  XCJ3  i 

COMMENT  6ET  Q  FOR  SHOCK > 

GA  *  -DELU  x  (CfiSG  »  ABStOELUi  +CONA  x  CSPCJ3)  /  VN  > 
IF  OA  <  0  AND  (CYCLECIO  OR  6tU.‘lCS3=i)  THEN  QA  0.01 
TLiMACJ3*OELX/<LIN£ARxCSPCJ3  ♦  C3SQ4  x  ASS (DELU. )  I 

IF  BURN! S3  J:  0  THEN  ELAST CPLAST XC  * 

IF  8URNCS3  s  3  THEN  BE6IN  EJ  *  EtJJ  I  VJ  *  VC J3  I 
•ATER  I  ENOI 

IF  BURN! S 3=4  THEN  BEGIN  VJ:=VCJ3»  PKASETRANSITION?  ENOI 
IF  BURHCS3  =5  THEN  STYROFOAM! 

IF  AtiS(PCJ3Kl.OO>3  THEN  PCJ>0I‘ 

VtJ3«VNI 

IF  TLIKACJXTLIMB  THEN'  BEGIN 
JCKIT*JI  TLIMB^TLIMACJJ  END! 

IF  IPCJ3  ♦  OC  J3)  >  PPEAK  THEN  BEGIN 

PPCAK  ♦  (PCJ3  ♦  OCJ3)  I  JPMAX  *•  J  I  END  I 
J*J*1>  J1*JH»  JWini 
IF  J*HESl3n  THEN  BEGIN 

IF  JXJSTAftn  THEN  GO  TO  GCTUVI 

COMMENT  TEST  VELOCITY  JUST  CALCULATED! 

IF  UCJSTARnJFO  THEN  JSTAR  *-JSTAR-*-l  END  ELSE  BEGIN 
IF  {CYCLE  MOO  X0)*0  THEN 
fcRITEtPOUTrFSUR’FACE)  ENOI 

IF  PLCELUJCC33SO  AND  JSTAR>CELU JCC3*8  THEN  BEGIN 
JPB  *  1  I  JPE  ♦  JSTAR  ♦  Z  I  RESTART  *  II 
I*  BOOL  THEN  GO  TO  FINISH  ELSE  GO  TO  READ  I  NO  J  '.-<01 
IF  T XMES2ABS (TWt NT 1)  THEN  BEGIN 
FLAP:  JPB  *  1  I  JPE  *  JSTAR  *  Z  l 

ENOI 

IF  (CYCLE  KOG  10)  SO  THEN  BEGIN 
KA1TE(P0MTCD6L3)I 
VRITE(POUTCOBL3)l 

*HIT£(POUT»F2.  TIN£(2>/«0»  JSTAR) I 

J*ll  MRITE(PCUT»FSVR»FACEIf 

J*  JPMAX  I  VRITE  (PO.'Tf  F3» ARRAYS)  * 

FOR  JC*I  STEP  1  UNTIL  JCC  00  BEGIN 
IF  J  >  JSTAR  THEN  GO  TO  KENT  I 
J*CELLtJC3»  MUTE  (POUT  »FGAOE»FOR*A6E»  END  I 


iAW f  -rt  ^  ’'vva/V  l,Am\,  *  ? 


TE  ' 


END  I 

NEWTl  0TNHX*0.6«TLIMS» 

IP  (DTNH1/DELT-1.1)>0  THEN  DTNHl*l.lxDELT» 

IP  DTNH1>DTMX  THEN  DTNH1*DTMX» 

OTN*DELTt  OEUT^OTNHl I  LIEUT I *DTN+DELT I  60  TO  LEPTI 

END  l 

pxnish:  end. 
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APPENDIX  n 

MATE  RIAL  DESIGNATIONS  AND  SUPPLIERS 


Designation  in  Text 
Aluminum  Powder*: 
Chunk 


Type  and-  Supplier 


MD801  Oversize  Aluminum  Powder. 
Metal*  Disintegrating  Corp. 
Elizabeth,  New  Jersey 

also 

Atoj.-'  zed  Grade  40  Aluminum 
Powder. 

Reynolds  Metals.  Co-. 

Louisville,  Kentucky 


Spherical 


Milled 


flake 


MD-AK  Aluminum  Foam 


Graphite  and  Carbon 


Silica  Microballoons 


No.  H-322  Aluminum  Powder 
Valley  Metallurgical  Co. 

Essex,  Connecticut. 

No.  A -547  Granular  Ball- Milled. 
Aluminum  Powder. 

Fisher  Scientific  Co. 

Fairlawn,  New  Jersey 

3xD  Aluminum  Powder 
Reynolds  Metals  Co. 

Louisville,  Kentucky 


MD-AK  Aluminum  Foam  .  . 

Emerson  and  Cuming,  XmL  ;  .... 
Canton,  Massachusetts  *. 

All  material  supplied  by  National 
Carbon  Company  and  referred  to 
in  text  by  manufacturer’s  designa¬ 
tions. 

Eccospfccres  SI 
Emerson  and  Cuming,  Inc. 

Canton,  Massachusetts  . 
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APPENDIX  m 

TABULATED  SUMMARY  OF  STATIC  AND  DYNAMIC  DATA 
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